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ABSTRACT 

The  propagation  of  small  distiirbances  in  an  infinitely  conducting, 
perfect,  compressible  mediiim  is  studied  from  the  standpoint  of  geometrical 
hydromagnetics.   The  mediiim  need  not  be  homogeneous.  A  wave-like  time- 
harmonic  distiirbance  is  assumed  given  on  an  arbitrary  reference  surface. 
Then  the  various  phase  fronts  (slow,  Alfven  and  fast)  which  may  be 
thought  to  evolve  out  of  this  stirface  are  obtained  as  solutions  first 
order  partial  differential  equations.   These  eikonal-like  equations 
are  solved,  as  in  optics,  by  means  of  rays,  which  are  shown  to  satisfy 
Hamilton's  equations  and  Fermat's  principle.   The  hydromagnetic  analogue 
of  Huygen's  phase  front  construction  is  described.   Next,  the  variation  of 
the  disturbance  strengths  along  the  rays  is  determined  in  terms  of  the 
disturbance  strengths  on  the  reference  siurface^to  complete  the  geometrical 
solution.   A  noteworthy  feature  of  this  solution  is  that  its  component 
modes  of  propagation,  slow,  Alfven  and  fast,  are  uncoupled  as  long  as  the 
medium  is  free  of  diffracting  objects  and  surfaces  of  discontinuity.   As 
an  illustration  of  the  general  theory,  the  complete  geometric  solution  for 
the  propagation  of  Alfven  disturbances  in  a  dipole  field  is  described  with 
the  aid  of  explicit  formulas.   Fast  wave  propagation  in  various  types  of 
inhomogeneous  magnetic  fields  is  also  briefly  discussed.   Finally,  the 
problem  of  reflection  and  refraction  of  slow,  Alfven  and  fast  disturbances 
at  (possibly  curved)  surfaces  of  discontinuity  is  treated.   The  hydro- 
magnetic  analogues  of  Snell's  laws  are  derived  and  their  relation  to  Huygen's 
wavelet  construction  is  elucidated.   It  is  shown  that  the  introduction  of 
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an  appropriate  set  of  direction  angles  for  the  unperturbed  magnetic  field 
(on  each  side  of  the  discontinuity)  l)  reveals  the  symmetries  of  "Snail's 
laws",  2)  leads  immediately  to  a  simplification  of  V. C.A.  Ferraro's 
laws  relating  to  the  reflection  and  refraction  of  Alfven  waves  and  3) 
makes  possible  a  simple  graphical  method  for  determining  the  angles  and 
speeds  at  which  the  multiply  reflected  and  refracted  waves  emerge. 


I.   IIWRODUCTION 

In  this  paper^  we  present  a  survey  of  propagation  phenomena  in 
magnetohydrodynamics  from  the  standpoint  of  "geometrical  hydromagnetics" 
The  medium  is  asstuned  to  be  a  classical  hydromagnetic  fluid  -  that  is, 
perfect,  compressible  and  infinitely  conducting.   The  imdisturbed  state 
of  this  medium  is  taken  to  be  time- independent;  however,  special  homo- 
geniety  is  not  required  so  that,  in  particular,  the  density  and  the 
embedded  magnetic  field  may  vary  from  point  to  point.   The  term  "geo- 
metrical hydromagnetics"  is  best  described  at  this  stage  by  saying 
that  it  bears  the  same  relation  to  the  linearized  equations  of 
magnet ohydrodynamlcs  as  the  geometrical  optics  of  crystals  bears 
to  Maxwell's  equations.   As  is  well-known,  there  are  essentially 
two  geometrical  theories.   One  involves  the  propagation  of  dis- 
continuities and  ptilses,  the  other  involves  time-harmonic  wave  motion. 
Since  the  propagation  of  discontinuities  has  already  been  treated  by 
J.  Bazer  and  0,  W.  Fleischman  [l]  we  confine  ourselves  here  to  discuss- 
ing in  detail  only  spacially  continuous,  time-harmonic  wave  motion. 
Then  at  the  end  of  Section  II,  we  indicate,  following  F.  G.  Priedlander 
[2a]  and  D.  Ludwig  [s]  how  to  include  both  theories  under  one  general 
formalism. 

Within  the  context  of  geometrical  hydromagnetics,  two  basic  prob- 
lems will  be  treated.   Suppose  that  the  disturbance  is  known  to  be  of 
the  form  exp  [icJW(x)-(t-t  )]  SV  (x,w),  -  «,  <  t  <  «^  on  a  fixed  reference 
surface  in  space.   In  this  expression,  ^   is  the  angular  frequency, 
5V(x,w)  is  a  real- valued  7-vector  whose  components  consist  of  the  dis- 
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turtances  in  density  and  in  the  components  of  the  magnetic  and  velocity 

fields,  W(x)  is  an  unknown  phase  fimction  andt  is  any  conveniently  chosen 
#•'  o 

reference  time.   It  is  to  be  understood  here  that  W(x)  =  0  is  the  equation 
of  the  reference  surface  at  time  t  .   Suppose,  in  addition,  that  w  is 
large  or  alternatively  that  the  wave  lenths  are  small  in  comparison  with 
the  significant  lengths  of  the  problem.   The  first  problem  is  to  deter- 
mine the  field  everywhere  when  the  undisturbed  medium  is  free  of  boundaries 
and  of  discontinuities.   The  second  basic  problem  involves  the  reflection 
and  refraction  of  these  small  disturbances  at  cirrved  surfaces  of  dis- 
continuity in  the  undisturbed  medium. 

Sections  II  and  III  are  devoted  to  the  solution  of  the  first  prob- 
lem. In  Section  II,  it  is  shown  how  to  decompose  the  basic  disturbance 
exp  [iwW(x)j6V(x,u))  on  the  reference  surface  into  three  pairs  of  funda- 
mental disttirbances  -  the  so-called  slow,  Alfven  and  fast  disturbances. 
Each  of  these  six  fundamental  disturbances  gives  rise  to  wave  motions 
of  the  form-^  exp  [iw(W^J^(x)  -  (t-t  )]5V^"^ '' (x,tj)  j   =  1,2,..., 6,  where  the 

W  "^  (x)  -  (t-t  )  =  0  are  the  equations  of  the  phase  fronts  which  originate 
^  o 

at  the  reference  front  at  time  t  =  t  .   The  anisotropic  character  of 

o 

wave  motion,  due  in  part  to  the  motion  of  the  undisturbed  medium  and  in 
part  to  the  presence  of  a  magnetic  field,  is  evidenced  by  the  fact  that 
the  propagation  speeds  of  these  phase  fronts  are  direct ion- dependent. 


We  neglect  terms  of  order  w"  and  higher.   It  is,  however,  possible 
in  principle  to  carry  out  the  analysis  to  arbitrary  order  in  w  . 


A  method  for  constructing  the  phase  functions  \J       ,    j  =  1,2,..., 6  is 
presented  in  Section  III  where  the  ordinary  differential  equations  for 
the  "rays"  associated  with  each  mode  of  propagation  are  derived.   The 
analogues  at  Huygens'  phase  front  construction,  Fermat ' s  principle, 
Presnel's  ray  surface  and  reciprocal  surface  of  phase  front  normals 

are  discussed.   As  in  the  geometrical  optics  of  crystals,  the  rays 

2 

are  the  conveyors  of  the  "energy  in  a  given  mode".   It  is  essentially 

this  fact,  specifically  the  conservation  of  this  "energy"  along  a 
tube  of  rays,  that  enables  one  to  determine  how  each  of  the  5V  "^  (x,  c^) 
behaves  as  x  varies,  thereby  completing  the  geometric  solution  of  the 
first  problem. 

The  fact  that  the  amplitudes  6  V{x,^)   can  be  determined  deserves  to 
be  underscored.   In  geometrical  optics,  for  example,  it  has  been  known 
for  a  long  time  that  the  transition  from  wave  optics  to  geometrical 
optics  involves  asymptotic  expansions  which  lead  to  methods  for  con- 
structing the  phase  fronts  such  as  have  been  described  above.  Not  so 
well  known,  however,  is  the  fact  that  the  geometrical  theory  enables 
one  to  determine  the  propagation  of  the  amplitudes.  A  systematic  ex- 
ploitation of  this  possibility  was  initiated  by  R.  K.  Limeberg  [U] 
only  twenty  years  ago  and  has  been  elaborated  and  extended  with  great 
success  in  the  papers  of  M.  Kline  [5a-c]  and  J.  B.  Keller  [6a-6c]  . 

As  in  geometrical  optics,  the  geometric  solution  can  be  obtained 
even  when  the  method  of  seperation  of  variables  fails.   Of  even  greater 
significance  in  practical  applications  is  that  each  of  the  six 
fundajuental  wave  motions  propagate  independently  of  one  another.  When, 
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by  contrast,  one  attempts  to  solve  linearized  equations  of  magneto- 
hydrodynamics  exactly,  one  finds  that  the  six  modes  will,  in  general, 
be  coupled.   Coupling  of  the  geometric  modes  does,  however,  take 
place  when  a  phase  front,  be  it  slow,  Alfven  or  fast,  is  incident  on 
a  surface  of  discontinuity  in  the  undisturbed  medium.   Then  one  finds, 
as  a  rule,  that  each  type  of  incident  wave  produces  three  refracted 
and  three  reflected  waves.   Section  V  is  devoted  to  the  discussion  of 
the  laws  governing  the  reflection  and  refraction  of  the  various  types 
of  incident  modes.   It  is  shown  that  the  introduction  of  an  appropriate  set 
of  direction  angles  for  the  magnetic  fields  (see  Figure  7)  makes 
possible  an  effective  graphical  procedure  for  determining  the  angles 
of  reflection  and  refraction  and  leads  to  a  simplification  of  V.C.A. 
Ferraro's  [7]  laws  of  reflection  and  refraction  for  Alfven  waves. 
¥e  also  explain  the  relation  of  this  graphical  procedure  to  Huygens' 
classical  construction  of  reflected  and  refracted  wave  fronts. 

The  theory  of  Sections  II  and  III  is  illustrated  in  Section  IV 
with  the  aid  of  several  concrete  examples.   In  all  cases,  the  medium 
is  assiimed  to  be  motionless  but  otherwise  arbitrary.   The  section 
begins  with  the  proof  of  somewhat  siorprising  property  of  Alfven  wave 
propagation,  namely,  that  the  strength  of  a  disturbance  on  an  Alfven 
front  is  independent  of  the  local  geometry  of  the  front.   Specifically, 
it  is  shown  that  the  strength  of  an  Alfven  disturbance  varies  as  the 
fourth  root  of  the  density  along  a  ray  and  is  independent  of  the  so- 
called  expansion  ratio  (see  equation  (3.21)).   This  result  means, 
among  other  things,  that  the  propagation  of  Alfven  fronts  into  a  region 
in  which  the  rays  converge  need  not  be  accompanied  by  an  intensification 


of  the  distrubance.   In  fact,  in  the  case  of  constant  density,  the  dis- 
turbance is  constant  along  a  ray.   Next  to  be  considered  is  the  problem 
of  Alfven  wave  propagation  in  a  dipole  field.   Here,  we  give  explicit 
formulas,  up  to  a  single  quadratirre,  for  the  motion  of  the  phase 
fronts.   This  quadrature,  it  may  be  added,  is  easily  evaluated  in  the 
case  of  constant  density.   Finally,  we  derive  formulas  for  the  pro- 
pagation of  fast  fronts  into  various  types  of  inhomogeneous  magnetic 
fields. 

By  introducing  suitable  "diffracted  rays"  and  systematically  apply- 
ing the  methods  of  J.  B.  Keller  and  co-workers  (see  reference  [6a] ),  one 
may  extend  geometrical  hydromagnetics  to  cover  the  effects  of  diffraction. 
Experience  with  the  corresponding  extension  of  geometrical  optics  (see 
[6a] )  indicates  this  theory  will  be  accurate  to  wavelengths  as  large  as 
the  diameter  of  the  diffracting  object.   The  application  of  the  geometrical 
theory  is  also  limited  at  the  high  frequency  end  of  the  spectrum,  owing  to 
the  failure  of  the  equations  of  hyiromagnetism.   Derivations  of  the  single- 
fluid  equations  with  a  scalar  conductivity  from  the  two-fluid  theory  point 
to  the  following  restrictions.   (l)  The  impressed  frequency  v  =  (u)/2jt) 
must  be  much  less  than  the  electron  gyro- frequency  v   ^   the  plasma  frequency 
V  ,   and  the  frequency  of  collision  v  between  the  two  species  of  particles j 
(2),  the  relation  v  <  <  v  must  hold.   Ihen  the  additional  requirement 
that  the  displacement  current  be  negligible  in  comparison  with  the  con- 
duction current  guarantees  the  validity  of  the  hyiromagnetic  equations. 
It  is  not  claimed  here  that  these  are  the  weakest  or  minimum  number  of 
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restrictions;  depending  on  the  mode  in  question  and  its  direction  of 
propagation  some  of  these  may  "be  relaxed  considerably  (see  references 
[26],  [27],  and  [28]). 

A  trief  survey  of  the  relevant  literature  is  given  at  the  end  of 
Section  II. 


-  7 


II.   THE  LINEARIZED  EQUATIONS  AND  THEIR  GEOMETRIC  SOLUTIONS 

In  the  absence  of  disturbances,  it  is  assumed  that  the  medium  is 
undergoing  steady  isentropic  flow  so  that  this  flow  is  governed  by 

V  •  H  =  0,  Aq 

V  X  (H  X  u)  =  0,  A^ 

p(u  •  V)u  +  a^  V  P  +  ^^  H  X  (v  X  H)  =  0,  A.    (2.1) 

V  •  (P  u)  =  0,  A 
S  =  const  (independent  of  x  and  t).  A, 

In  these  equations,  H  is  the  magnetic  intensity,  u  is  the  fluid  velocity, 

P  is  the  density,  a  is  the  speed  of  sound,  V-  is  the  specific  magnetic 

inductive  capacity,  and  S  is  the  entropy  per  -unit  mass.   MKS  imits  are 

used  throughout.   We  prefer  to  ignore  gravitational  terms  here  and  in 

the  body  of  our  analysis.   However,  in  Section  IIZB  we  list  the  few 

chajiges  in  the  final  results  that  are  necessitated  by  the  inclusion  of 

such  terms.   The  ( time- dependent )  linearized  system  is: 

0  =  Lq(5H)   =  V  •  5H,  5Aq 

0  =  L,  (5H,  6u)  =  5R  +  y  X  (5H  X  u)  +  V  X  (H  X  6u),       6A., 

0  =     L^  (S  P,  5u,  SH)   =  P  5u^  +  P(u  *  v)5u  +  P  (5u  •  y)u   &A^ 

+  H  H  X  (y  X  5H)  +  V  (a^  SP)  (2.2) 

+  5P(u  •  V)u  +  U5H  X  (v  X  H), 

0  =  L^  (6P,  6u)  =  5P,  +  V  •  (u  SP)  +  V  •  (P  5  u),         5A 
3      •^      t       -  -  i 

0  =  &S  .  5A^ 

In  these  equations,  5H,  Su,  Sp  and  5S  denote  small  variations  of  H,  u,  p  and 

S.  L„,  L- ,  L^  and  L^  are  abbreviations  for  the  linear  differential  operators 
0'  »-l'  «»^2      3 


appearing  in  the  last  members  of  5A^  SA.^  ^  SA^  and  5A  ;  the  subscript 
t  signifies  partial  differention  with  respect  to  time. 

Forming  the  relation 

2 
^J6H  •  L^  +  Su  •  L^  +1  SP  L^   =  0,  (2.3) 

"by  multiplying  bA  ^bA     and  5A     by  6H«,    &u*  and  a  5P/P^   we  obtain  a^ter 
some  manipulation  the   "energy"  equation 

0     =    C    +  V  •    F  +  P  5  u  •  (5u   .   V)u  -  -     SP  6u   •    Fh  X  (v  X  H)l 

2 
+  H6u  •  r6H  X   (v  X  H)l    -  i  ~  u   .   V  (a^)    +  J   HoH^V  u  (2.1l.) 

-  US  H   •     (5H  •    v)  u, 
where  the   "energy"      C   is  defined  by 


£ 


2 
=  i  P  5  u^   +  i  ^i6H^+  i  |-5p2,  (2.5) 


and  the  "flux  of  energy"  or  system  "flux"  is  defined  by 

F  =  Cu  +  a  5P  6u  +  ^^  (H  X  5u)  X  5H.  (2.5') 

The  last  term  in  equation  (2,6)  will  be  recognized  as  the  Poynting 
vector,  since  in  fluid  of  infinite  conductivity  the  electric  field  is 
given  by  the  expression  M-HX  5u.   The  first  term  in  the  right  member  of 
(2.6)  gives  the  energy  density  transported  due  to  the  motion  of  the 
fluid  and  the  second  term  is  the  ordinary  flux  term  of  acoustics. 

Geometrical  hydromagnetics  enters  the  pictirre  when  one  asks  for 

a  solution  of  the  system  (2.2)  which  is  a  superposition  of  waves 

jij  — - — -^ 

C  should  not  be  confused  with  the  second  order  energy  term  obtained 

by  expanding  the  expression  for  the  total  energy  of  the  nonlinear  system. 

In  addition  to  terms  such  as  those  appearing  in  C  ,  the  second  order 

energy  contains  terms  involving  second  variations. 
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of  the  form 

00 

(5H,  Su,  5p)  =  e^^t^^i^H-tJ]  X  i^)''   (6H^"\  &u^"^  6p )      (2.6) 

n=0 

Here  w  denotes  the  angular  frequency,  W(x)  denotes  an  'unknown  phase 

function  and  6H^  ,   5u^  ^>  &P   ^  n  =  1,  2,  . , . ,  denote  unknown 

amplitude s .   These  amplitudes  are  required  to  be  real  apart  from  a 

common  phase  factor  of  unit  amplitude  and  to  "be  independent  of  w  and 

t.   It  is  not  required  that  the  infinite  sum  converge;,  merely  that  it 

be  asymptotic  in  some  sense  to  the  exact  solution  for  large  w  or  small  ^. 

For  a  fixed  instant  of  time, 

W(x)  -  (t  -  t  )   =  k 
-■         o       o 

is  for  each  constant  k  the  equation  of  a  surface  of  constant  phase  or 

phase  front.   To  complete  the  formulation  of  the  problem,  it  is 

.„  c--rr(n)   ~  (n)    -.   ^An)  _   ^  ^ 

necessary  to  specify  oH   ,  qji    and  oP    n  =  0,  1,2,  ...,  on  some 

conveniently  chosen  phase  front.   We  choose  for  this  purpose  the  phase 

front  A{t    )   defined  by 


J 


(t  )  :   W(x)  =  0  (2.7) 

o       ^« 


and  write 

n     (n)    (n) 


(6H,  5u,  5P)       =  2_   (^)   (5H   ,  Su   ,  5p  ^^0 

♦    "  L  r\        1^        O      "-O   '     O 


(2.8) 
n=0 
(to- 

.4(t  J. 


for  the  form  of  the  disturbance  on  .„, 

o 


One  may  approach  the  problem  of  determining  the  unknown  amplitudes 
and  phase  function  W(x)  in  either  of  two  ways.  One  may  seek  to  deter- 
mine  the  unknown  amplitudes  on  fixed  phase  fronts  W(x)  =  const  in  space 
or  follow  a  given  phase  front,  e.g.  yg(t)  defined  by 

(t):   W(x)  -  (t  -  t^)  =  0,  (2.9) 
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from  the  time  t  that  it  evolves  out  of  B«^(t  )  and  seek  to  determine 
the  unknown  amplitudes  on  «^(t).  We  shall  adopt  the  second  point  of 
view. 

As  already  observed  in  the  introductory  section,  an  arbitrarily 


polarized  wave  on  ^(t  )  gives  rise  to  three  pairs  of  phase  fronts,  slow, 
Alfven  and  fast  .  Each  of  these  six  fronts  carries  a  disturbance  of  the 
form  (2.6).   Since  the  amplitudes  of  each  distrubance  wave  can  be  deter- 
mined independently  of  the  others,  we  automatically  achieve  the  separa- 
tion of  wave  motions  referred  to  in  Section  I. 

The  velocity  of  the  phase  front  ■~^jit)   is  defined  to  be  the  velocity 
along  the  orthogonal  trajectories  to  the  family  of  phase  fronts, 
W(x)  -  (t-t  )  =  0.   Employing  this  definition  one  finds  that 

U  =  n/  |V  W]  (2.10) 

where  n  is  the  unit  normal  to  the  phase  front  so  that 

n  =  W  /  |VW|.  (2.11) 

If  we  substitute  the  right  member  of  equation  (2.7)  into  equations 
(2.2)  BA_  -  (2.2)  6A  and  equate  the  powers  of  w,  we  obtain  the  follow- 


ing recursive  system: 

3h(0)  =  0, 

(2.12) 


n  .  5h(°^  =  0, 


(M  -P^  c  I)  5R^°^   =  0; 

n  •  5H^^^   =  -  UL  m^^-^h, 


1 


,(k)  _  x>/^o(k-l) 


(2.13) 


(M  -p  2  c  I)BR^''^   =  P(6R^'^"^0,   k  =  1,2,,  . 
In  equation  (2.I3),  M  is  the  symmetric  6  >^  6   matrix 

Let  X  =  X  (t)be  one  such  trajectory.   Then  we  find  that  U  =  ;r*'  and 

W(x(t))  -(t  -  t  )=  0  (identically).   Differentiating  this  last  equation, 

-1 
we  conclude  that  U  =  |VW| 


n 


M 


0 

y 

n 

0 

0 

0 

1 

2 

H  H 

0 

0 

aP* 

0 

0 

y 

!_ 

M-^H 

0 

0 

0 

0 

0 

n 

1 

0 

ap2 

0 

0 

0 

0 

~ 

—      —  — 

—     —  _ 

J 

—  - 

~~l_ 

0 

0 

0 

0 

0 

-I^^H 

0 


0 


0  I  -M-^H 


(2.14) 


5R         is  the  colxiran  vector 
,(k) 


5R' 


p   y     n'   yP       z    z 


;  (2.15) 


P(5R     )  is  the  column  vector 


P(6R^^"^^) 


r   i_(k-l)     _(k-l)  _i_(k-l)  ^  _(k-l) 

i  (k-i) 


(2.16) 


1  _(k-l)     U   2 

"^'       (^)   L- 


2.Z 


'P 


'l,z 


and 


c   =  U  -  u 


n 


(2.17) 


is  the  disturbance  speed  (relative  to  the  flow  velocity  in  the  direction 


normal  to  the  wave  front).   The  operators  L  , 


L  appearing  in 


equation  (2.l6)  are  the  same  as  the  L^  ...  L  that  appear  in  (2.2)  ex- 
cept that  the  time-derivative  operators  have  been  deleted.   The  second 
subscript  on  the  vector  operators  indicates  which  component  is  intended. 

The  particularly  simple  form  of  M  is  the  result  of  choosing  a 
local  system  of  coordinates  on  (j/(t)  such  that  the  x-axis  is  directed 
along  n  and  the  y-axis  is  directed  along  H  -  (H  •  n)n. 
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We  shall  henceforth  he  concerned  only  with  the  lowest  order  terms 
in  the  solution  (2.6).   This  means  that  we  are  required  to  find  the 
solutions  of  the  zero-th  order  equations  (see  equation  (2.12)).   With 
the  aid  of  the  remaining  equations  one  may,  in  principle,  obtain  the 
first  and  higher  order  terms;  however  the  problems  that  one  encounters 
in  practice  do  not  as  yet  seem  to  warrant  such  refinement. 

We  begin  by  observing  that  5R    is  necessarily  an  eigenvector 

of  M.   Since  M  is  symmetric  it  has  a  complete  set  of  six  real-valued 

mutually  orthogonal  eigenvectors.   To  find  these,  we  first  determine 

i_ 
the  eigenvalues  P^c  from  the  seciilar  equation 

0  =  det  (M-P^cl)  =   (Pc^-H^)  [(Pc^-  Pa^)(Pc^-^lHn^)-Pc^(H2-H^)]  .  (2.18) 
The  solutions  of  this  equation  evidently  occur  in  pairs  +  c,  c  >  0. 
Let 


c.   =  b 
A      n 


(^^  H  ^  P"^)' 
n 


(2.19) 


-  the  so-called  Alfven  disturbance  speed  -  be  the  non-negative  zero 
of  the  first  factor   and  c  ^{Oslow)   and  c„(c^j^g^)  be  the  smallest  and 
largest  of  the  non-negative  zeroes  of  the  second  factor,  in  the  last 
member  of  equation  (2.18).   Then,  it  is  found  that 


I  (1  +  r)  -  i  [(1  +  r)  -  Ur  cos  x] 


i\ 


=_i  _ 


I  (1  +  r)  +  i  [(1  +  rf   -  J+r  cos^X]' 


2 


(2.20) 


(2.21) 


and  that 


0  <  Cg  <  min|  a,  bj^Z  <  max  |  a,  b^l  <  c^, 

where  r  and  b  are  defined  by 
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r     =     a  /b 


^     =  i/l^  H^  p-1 


(2.22) 
(2.23) 


and  X.  is  the  angle  between  H  and  n.      Figure  1  is  a  polar  diagram  of 

w 

Gg/b,  cj^   and  c  /b  against  \   The  figure,  which  is  plotted  for  the 
parameter  vulue  r  =  l/2  is  typical  of  those  in  the  entire  parameter 
range  r  <  l/2.   Representative  curves  associated  with  values  of  r  in 
the  ranges  l/2  <i-<l,  r  =  l,  l<r<2,  2<r  are  given  in  the 
report  [l  aj  . 

Let  it  be  supposed,  for  the  present,  that  n  is  neither  parallel 
nor  perpendicular  to  H.   Substituting  the  appropriate  value  of  c  into 
equation  (2.12)  and  solving  for  the  components  of  5R   ,  we  find  that 
6H   ,  5u    and  SP    in  Alfven,  slow  and  fast  waves  are  given  by 
the  following  formulas: 

Alfven: 

1 
-,2 


m[^K    +  (e/E^)  n  X  H, 


Su 


(0) 


e((i  H  /E^Pb  )  n  X  H, 
n'    n  --   .-v' 


(2.2^) 


6P^°^  =  &S^°^  =  0; 


E 


J~P    b  Isin  X 


Slow  and  FsLst : 


6H^'^-' 

=   -    (e/E^ 

)n  X    (n  X  H), 

_    e 
=   +    — 1 
E  2 

r          1.2 

r    r  n        ,  V 

g(      2   ~   ^'^ 
c 

5    (0) 
Op 

=  -     4    P 
E^ 

e 

PC 


,  (2.25) 
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E^  =  /p  b|sin  X| 


2a  2 


2(c7T3^)-(14  r) 
Cc2/b2)-  r 

2^  J  |(c%2)  -  cos^Xj 


c 


(c%2)_(i+^)) 


In  equation  (2.25),  c  is  c  or  c  according  as  to  the  wave  is  slow  or 
fast.   There  are  two  solutions  associated  with  each  of  the  speeds  c., 

c  and  c„;  one  corresponding  to  motion  of  the  phase  front  along  n  (upper 

s      f  ^ 

signs)  and  the  other  corresponding  to  motion  in  the  opposite  direction 
(lower  signs)  along  -n.   The  E's  are  normalization  constants,  having  the 
dimensions  of  energy  per  unit  volume,  which  are  fixed  by  requiring  that 
p/2  times  the  square  of  the  Euclidean  length  of  the  vector  5R    be  equal 
to  e  -  in  symbols  (see  equations  (2.15)  and  (2.5)) 

ip^^iO)    .  .^(0)  ^  i(p  ^J0)^2  ^  i  ^  (5^(0)^2  ^  i  ^  ^^^(0)^2  ^  ^2^  (2.26) 

The  quantity  e,  which  is  evidently  dimensionless,  is  called  the  strength 

2 
of  the  wave.   Measured  is  units  of  E,  €  is  the  energy  per  unit  volume  in 

the  wave  (see  footnote  3).   Ii^  order  that  the  linearization  on  which  the 
entire  analysis  is  based  be  valid,  it  is  necessary  that  |e|  be  a  small 
quantity  initially  and  remain  so  subsequently. 

The  analogues  of  equations  (2.2^+)  and  (2.25)  foi"  the  thus  far  ex- 
cluded directions  of  n  may  be  obtained  as  follows:  Rotate  n  into  the 
direction  parallel  (or  normal )  to  H  keeping  the  direction  of  n  X  H 
fixed.  Then  calculate  the  limits  approached  by  the  expressions  for 
5H^   ,  5u^  ^   and  5p ^   .   The  resulting  limits  are  the  desired  formulas 
as  may  be  verified  by  substituting  the  appropriate  components  into 
the  eigenvalue  equation  (2.12). 
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By  direct  calculation^  it  can  be  shown  with  the  aid  of  equation 
(2.18)  that 

i  p(5u^°^)2  =  i  4  ^'^^^°^)'  -^  -   ^   (^H^°^^'  =  *  ^'  ^2.27) 

in  all  waves,  slow,  Alfven  and  fast.   Thus,  the  energy  density 
associated  with  each  of  the  fundamental  modes  is  divided  in  equal  parts 

"between  the  kinetic  energy  density  on  one  hand  and  the  magnetic  and  in- 

5 
ternal  energy  densities  on  the  other  . 

Suppose  n  X  H  :[=  0  and  that  n  and  H  are  in  the  plane  of  the  page. 
Then,  a  glance  at  equations  (2.2U)  and  (2.25)  reveals,  the  magnetic 
vector  5H  is  polarized  perpendicular  or  parallel  to  the  plane  of  the 
page  according  as  the  front  is  an  Alfven  front  or  not.   In  all  cases, 
fJH  is  parallel  to  the  phase  front  and  its  direction  is  determined  by 
n  and  H. 

The  problem  of  determining  the  disturbance  evolving  out  of  a 
given  initial  disturbance  may  now  be  reduced  to  the  corresponding 
problem  for  the  six  fundamental  modes.   Let  x  be  a  point  on  the  phase 
front  J  it    )   and  let  n(x)  be  the  unit  normal  at  x.   Corresponding  to 

each  such  x  and  n(x),  let  R^  ^  =  R^  q^^^'Q) ,    ^^,0  "  ^A,0  ^^^2.^'  ^f  " 

+ 
R   (x,n)  be  the  six  mutually  orthogonal  eigenvectors  of  the  matrix 
f ,  0  •  ** 

M.  Without  loss  of  generality  it  may  be  assumed  that  these  vectors 
are  of  unit  length  .  Then,  ignoring  all  but  the  lowest  order  terms 
in  the  right  member  of  (2.8),  we  may  decompose  the  initial  disturbance 

■^The  internal  energy  density  vanishes  in  Alfven  waves. 

The  components  of  these  unit  vectors  may  be  obtained  from  equations 
(2.2ii-)  and  (2.25)  simply  by  setting  6=1. 


16  - 


into  the  six  fundamental  waves  as  follows 

aSP 


5D(x) 


1 

M.2 


Q     6H  ,  Su  ,  Su  ,  ^    Su  ,  (^)'  6H 
P     y    n    y   P     z   P     z 


E  (5D't(x)  +  5D".(x)) 


(2.28) 


J(- 


Here  x  is  any  point  on  ^ (t  )  and  j  takes  on  the  "values"  s,  A  and  f  and 
5irT(x) 


+    +  +  + 

er  „  R7  ^.   Since  the  RT  „  are  unit  vectors,  the  £-   _  measure 

the  strengths  of  the  fundamental  waves  in  the  ahove  decomposition.   If 

6D(x)  is  complex- valued  then,  as  the  Rr  ^  are  real,  the  e?  ^  will  he 

complex-valued.   However,  we  shall  show  at  the  end  of  Appendix  I  that 

only  the  amplitudes  of  the  er  ^  change  on  the  phase  fronts  -,(^":(t) 


j  =  s.  A,  f,  which  evolve  out  of 


^\(t  J. 


The  distinction  between  real 


and  complex  e-  is  therefore  not  an  essential  one. 
3 

We  have  thus  far  ignored  what  may  be  termed  contact  fronts.   Con- 
tact fronts  do  not  propagate  relative  to  the  londerlying  flow.   In  addition, 
they  may  carry  a  disturbance  6S  in  the  entropy  S.   Wo  special  consideration 
is  required  for  contact  waves  that  appear  as  limiting  cases  of  the  slow, 
fast  and  Alfven  modes.   However,  if  it  is  desired  to  take  into  account 
contact  fronts  that  are  not  connected  to  the  fundamental  modes  by  means  of 
of  a  limiting  process,  it  is  necessary  to  increase  the  number  of  com- 
ponents in  6D  and  RT,  3   =     s,  A,  f  by  one  and  to  add  an  entropy  vector 
e  R  to  the  sum  (3.28).   To  obtain  the  components  of  R  one  simply  sub- 
stitutes the  right  member  of  equation  (2.8)  into  (2.2)  and  looks  for  a 


solution  in  which  c  =  U  -  u  vanishes. 

n 
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Much  of  the  work  on  the  propagation  of  continuous  time -harmonic 
disturbances  in  a  compressible  hydromagnetic  medium  has  been  con- 
cerned with  the  derivation  of  exact  solutions  for  various  simple 
geometries.   The  contributions  of  N.  Herlofson  [S^ ,  H.  C.  van  de  Hulst 
[9],    5.    Lundquist  [lO]  ,  A.  Banos,  Jr.  [ll]  ,  to  mention  a  few  of  the 
earlier  papers,  are  examples  of  such  work.   Herlofson  and  van  de  Hulst 
were  the  first  to  note  the  existence  of  essentially  three  distinct 
disturbance  speeds  (of  equations  (2.I9)  -  (2.21)).   The  more  recent 
asymptotic  results  of  M.  J.  Lighthill  [l2]  are  geometrical  in  character 
but  valid  only  when  the  undisturbed  state  of  the  medium  is  spacially 
homogeneous  as  well  as  time-independent.   In  the  reports  [13]^  [?-''il  > 
K.  0.  Friedrichs  and  H.  Kranzer  called  attention  to  hyperbolicity  of 
the  linearized  hydromagnetic  equations  -  a  property  that  it  shares 
with  the  linearized  equations  of  gas  dynamics  and  Maxwell's  equations. 
In  addition,  they  took  the  first  steps  in  exploiting  the  analogies 
between  these  systems.   Among  other  results  they  described  the  wave 
forms  emanating  from  a  weak  localized  disturbance  in  a  uniform  medium 
and  catalogued  the  various  fundamental  wave  motions.   The  propagation 
of  weak  hydromagnetic  disturbances  was  discussed  by  J.  Bazer  and  0.  ¥. 
Pleischman  in  their  report  [la]  and  paper  [lb] .   It  should  be  noted 
that  this  paper  dealt  with  propagation  in  inhomogeneous  media.   In  [15] 
H.  Grad  discusses  continuous  wave  motion  in  a  uniform  media.   Among 
other  things,  he  shows  that  in  more  than  one  dimension  the  fundamental 
waves  can  only  be  partially  separated,  in  contrast  to  the  complete 
separability  of  the  geometrical  modes.   Grad  shows  however,  that 
special  solutions  can  be  constructed  which  propagate  at  one  of  the 
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above  mentioned  disturbance  speeds  (relative  to  the  flow).   The  reader 
is  referred  to  the  review  article  by  G.  F.  J.  MacDonald  [l6i]  for  a 
survey  of  geophysical  applications  of  hydroraagnetic  wave  propagation 
and  for  additional  bibliographical  material.   The  paper  [17]  by 
A.  J.  Dessler  and  E.  W.  Parker  on  the  hydromagnetic  theory  of  geo- 
magnetic storms  should  also  be  consulted  in  this  connection. 

A  comparison  of  the  formulas  (2.2^)  and  (2.25)  for  the  six 
fundamental  propagating  modes  with  the  corresponding  formulas  for  the 
jumps  across  weak  shocks  (cf .  [lb]  p.  268  and  [13]  pp.  1^1-20.  )  reveals 
that  they  are  the  same  apart  from  a  constant  multiplier.   Tliat  this  is 
to  be  expected  may  be  seen  from  the  following.   Note  that  the  f\mda- 
mental  asymptotic  expansion  (2.7)  has  the  form 

(5H,  &u,  Z^)      =  E    f^(cp  (x,t))  (5H^"^,  Bu^''^  6P^^^)       (2.29) 
n=0 

where  ^(x.t)   =  W(x)  -  (t  -  t  )  and  (2.3O) 


f   =  exp  (iws) 

f   =   /  f  ,(s)  ds 
n     ./    n-1 


(2.31) 


Wow  it  is  easily  verified  that  the  right  member  of  equation  (2.29)  is  a 
formal  solution  of  equation  (2.2)  for  arbitrary  f  (s);  provided  f  is 
defined  as  in  (2.3I).   This  function  may  in  fact  be  highly  singular^  e.g.  a 

delta  function  or,  more  generally,  a  distribution  in  the  sense  of 

[29] 
L.  Schwartz    .If^in  particular,  one  chooses  f  (s)  to  be  a  function  which 

takes  the  value  1  within  the  boundary  of  a  disturbance  and  zero  outside, 

one  is  led  to  the  theory  of  Bazer  and  Fleishman  [la,  b]  and  to  the 

formulas  (2.2^)-)  and  (2.25)  for  the  fundamental  propagating  modes. 
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Formal  expansions  of  the  above  type  were  first  introduced  by 
F.  G.  Friedlander  in  his  study  [2a]  of  propagation  phenomena  in 
acoustics.   D.  Ludwig  [3]  has  employed  similar  expansions  in  his 
treatment  of  the  initial  value  problem  for  hyperbolic  systems.   This 
work  contains  the  most  comprehensive  treatment  to  date  of  various 
mathematical  questions  which  arise  in  the  rigorous  theory  of  such 
expansions.   The  earlier  papers  by  R.  Courant  and  P.  D.  Lax  [I8] 
R.  Lewis  [19]  and  P.  D.  La>:  [20J  ^  (particularly  the  last  of  two)  also 
bear  directly  on  the  problem  under  discussion. 

Much  of  the  material  described  in  the  last  two  paragraphs  may  be 
regarded,  as  an  outgrowth  of  the  work  on  geometrical  optics  by  R.  K. 
Luneberg  l^] ,    on  geometrical  acoustics  by  J.  B.  Keller  [6bJ  and  on 
general  hyperbolic  equations  by  M.  Kline  [5j.  The  present  work  is 
to  a  large  extent  patterned  after  the  work  of  these  three  authors. 
An  informative  historical  account  of  the  evolution  of  the  geometric 
method  has  been  given  by  M.  Kline  [5c] .   Among  other  things,  the 
relation  of  the  geometrical  method  to  the  W.K.B.  method  is  explained 
in  this  paper. 

Finally,  it  should  be  mentioned  that  M.  J.  Lighthill  [l2]  ,  F.  G. 
Friedlander  [2b]  and  H.  Weitzner  [2l]  have  derived  expressions  for 
the  far  fields  due  to  various  types  of  point  sources  in  a  homogeneous 
medium. 
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III.   THE  GEOMETRICAL  SOLUTION  COMPLETED 

A.   Rays  and  Phase  Fronts 
Setting 

p  =  VW,        p  =   |V  W|  (3.1) 

we  find,  employing  the  definitions  of  n,  U  and  c  (see  equations 
(2.10),  (2.11)  and  (2.1?)),  that 

n  =  p"^p  ,      U  =  p"-*-  (3.2) 

c  =  p"  -  (u  •  p)p'  (3.3) 

The  vector  p  is  called  the  phase  normal. 

In  this  notation,  we  may  re-express  the  determinantal  equation 

(3.3)  as 


0  =  }f(x,p)  =  ri  >li  (x,P,)  ?Ji(x,p),  (3.M 


where 

+ 


The 


T=i  .   (x,p)  =   (u  •  p  -  1)  +  p  c.  (x,  n). 

1=1,2,3.  (3.5) 

quantities  c^ ,  c^  and  c_  are  the  disturhances  speeds  c.,  c  and 
c  respectively.   The  functions  p,   /f-,-,    n   2   ^^^^^  ^  will  he  referred 
to  as  the  total,  Alfven,  slow  and  fast  Hamiltonians  respectively  for 
reasons  which  will  be  clear  presently.   Since  by  definition  p  =  V  W, 
each  of  the  six  equations 

^i(x,p)  =  0  (3.6) 

is  a  first  order  partial  differential  equation.   As  is  well  known,  the 
problem  of  solving  such  equations  may  be  reduced  to  that  of  solving  a 
system  of  ordinary  differential  equations  of  the  first  order  in  (x,p) 
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space.   The  rays  are  the  curves  in  x- space  obtained  hy  solving  this 
system  of  equations.   If  the  time  t  is  employed  as  the  curve  para- 
meter then  the  ray  equations  corresponding  to  each  of  the  above 
Hamiltonians  are  simply  Hainiltons'  equations  of  motion  with  x  and  p 
serving  as  space  and  momentum  variables^  respectively.   Employing 
this  fact  we  are  led  to  the  following  set  of  ray  equations  (cf 
[lb]  p.  371): 
+  Alfven  rays : 

(a) 


J.  =  u  +  [\^/Pr   sgn  (H  •  p)  H  =v  -, 
dt    "  "  -   '"  '^ 


d| 
dt 


(p  •  V) 


X  Jt  X 


+  slow  rays: 


dx 

dt 

dp 
dt 


u  +  (m./p)^  sgn  (H  •  p)  H 

[u  +  (H/Pf  sgn  (H  •  p)  h]   , 

■J 

-  -  s  u.  -  c^{y.,n)Qi[Y.,n)    -         n'-'  --s 

=  -(p  •  v)  u  -  p  X  (\7  X  u) 


(b) 


+P  Vc^, 


(a) 


(b) 


+  Fast  rays 
dx 


=  u  +  C4:>  n  + 


dt 


_   b^r  (Hth/H^)     ,„   „   ,    t 

"^  — 7 \   o.t     — T   (H  -  H  nj  =  v^ 

(x,n)  Q(x,n)   ^-    n      -f 


^  _ 


-  (p  •  V)  -  p  X  (\7  X  u)  +  p  V  c^ 


(a) 


(b) 


(3.7) 


(3.8) 


(3.9) 


The  upper  and  lower  signs  in  the  right  members  of  these  equations  cor- 
respond to  the  "+"  and  ''-"  Hamiltonians,  respectively.   In  equation 
(3.7)  sgn  (h  •  p)  denotes  the  sign  uf  H  •  p  and  in  equations  (3.8)  and 
(3.9)  the  quantity  Q  is  defined  by 


-»*»    yty 
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{1  +  rf  -   k  r    (R  •    nf  /   H^ 


2   , 


(3.10) 


the  ratio  is  defined  in  equation  (2.22).   The  vectors  v.   (v.  ). 

•^  A    —A 

V  ,  (v  )  and  Vj,  (v„  )  are  called  the  forward  (haclcward)  ray 
~s^   »-s       *•  f  Vf  '  -^- —  — — 

velocities.  Observe  that  in  all  cases 

dx  ,    , 

p  •  —  =  1;       t  >  t  (3.11) 

~   o 
"-   dt 

where  t  is  regarded  as  the  initial  instant.   This  relation  follows 
o 

directly  from  the  ray  equations  and  the  fact  that  u^p  +  p  c.  =1, 

i  =  1,2;,3  [see  equations  (3.5)  and  (3.6)3-   It  should  be  observed  that, 

as  in  crystal  optics,  the  phase  normal  p  is   not  in  general  directed 

along  the  ray. 

Now  let  X  =  f   (2,  ;^  )  be  the  equation^,  in  parametric  form,  of 

any  member  of  the  family  of  phase  fronts  (see  equation  (2.9  )  and  let 

p  =11  (I  ,^  )  the  wave  normal  at  x  =  f  (|^,|^).   Then  one  can  in 
w-^12  v-'^J-^ 

principle  determine  solutions  of  each  of  the  above  ray  systems  of 
the  form 

^i-   =   Ji-  (^r^2^  *  -  ^o)  (3.13) 

such  that  f.-  (^-L^lg'  °^  =£°(^1^^2^  andJI.-  H^,^^,    0)  =  n° H^,^^) 


'Note  that  p°  =  p°n,  where  n  =  i:i(^-,  ,|p)  is  the  unit  normal  to  the  phase 
front  at  i^^,^^).      Substituting  x  =  I!  {^j_>^2^  ^^^  ^     =  p  n  into 

n.-   (x,p)  =  0,  i  =  1,2,3  we  obtain,  in  general,  six  values  of  p  there- 
by  fixing  the  initial  value  of  p  in  each  of  the  equations  (^4^.12)  -  (i+.li+). 
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for  all  i  .   It  is  assumed  here  and  in  equation  (3.12),  that  1  takes  on 
the  values  s,  A  and  f.   Equations  (3.12)  and  (3.8)  expresses  in  para- 
metric form  the  six  phase  fronts  evolving  out  of  the  initial  manifold. 
In  geometrical  terms  the  situation  maybe  described  as  follows:  6j(t  ), 
in  Figure  2,    represents  a  small  neighborhood  of  the  point  x°  =  f°(P      P    ) 

(t^).   This  elemental  phase  front  splits 
up  into  six  elemental  phase  fronts  6  J.-(t),  i  =  s,  A,  f,  propagating 
with  the  instantaneous  velocity  y. -,  i  =  s.  A,  f  along  their  associated 
rays.   The  positions  and  orientations  of  the  6  J.i(t)  at  each  instant 
of  time  is  given  by  equations  (3. 12)  and  (3. 13),  respectively.   As  the 
parameters  |^  and  ^^  range  over  their  domain  of  definition  5j^(t  ) 

sweeps  out  the  reference  phase  front  Ji{t    )  while  the  bj-(t    )   sweev 

0  °  i  '  o'     -^ 

out  the  six  fronts  d^-(t),  into  which  J(t  )  has  evolved. 
+ 
^et  /i~  ^  (i  =  s,  A  and  f)  denote  the  flux  (see  equation  (2.6)) 

associated  with  a  fundamental  wave  of  type  (+,  i).   Combining  equation 
(2.6)  with  the  formulas  for  5H^°\  bj^\    5P^°^  and  simplifying  the 
results  with  the  aid  of  equation  (2.I8),  it  is  found  after  some  manipula- 
tion that 

■ni  +     r      +\2    + 

F.-  =  (e.-)  V.-,        i  =  s,   A,   f  .        (3.14) 

In  this  equation  ;,  (e^-)   is  the  energy  density,  measured  in  units  of 
E^,  (see  last  equations  of  (2.24)  and  (2.25))  on  the  (+,  i)  -  ray  at 
time  t.  An  explanation  of  how  to  calculate  the  e  -  from  their  initial 


Solving  for  l^,!^  and  t-t^  in  terms  of  x^,  x^,  x  in  equation  (3. 12)  we 
obtain  three  equations  one  of  which  is  of  the  form  t-t  =  W(x).   In  this  way 
way  we  obtain  W(x). 
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We  shall  henceforth  append  the  subscript  i,  i  =  s,A,f  to  the  quantities  F,e, 
E,  etc.  in  order  to  identify  the  fundamental  wave  to  which  these  quantities 
correspond. 


2k 


values  (see  equation  (2.28))  will  be  given  in  the  next  section.   Here, 
ve  merely  note  that  equation  (3.1^)  is  the  precise  form  of  the  state- 
ment "the  rays  are  conveyors  of  energy". 

B.   Variation  of  the  Strengths  of  the  Fundamental 
Waves  along  the  Rays. 
Let  ^.—,   t  >  t  denote  the  (+,  i)-th  phase  front  that  has  evolved 
out  of  the  reference  phase  front  ^(t  )  (see  equation  (2.7)).   Let 

5  D  i(x)  =  €-(x)  R,i{x,n)  (3.I5) 

1  "^      1  -   1  ~  •> 

denote  the  disturbance  at  a  point  x  of  >C'-~(t).   As  in  Section  II, 
n  =  n(x)  is  the  normal  to  ^ .-{t)   at  x  and  R.-(x,n)  is  the  normalized 
eigenvector  of  the  matrix  M  (equation  (2.1^+))  associated  with  the  dis- 
turbance speed  +  c..   To  complete  the  geometrical  solution,  it  is 
necessary  to  determine  the  -variation  of  6D.— (x)  at  each  point  of 
«/}.  -  (t )  for  all  t  >  t  ,  given  that  5D-(x)  is  specified  on  jj'(t  )  by 

=  4  M")   R?  .(?^n),  (3.16) 


+, 


&D.-(x) 


^     1,0     1,0 
o 


(see  equation  (2.28)). 

Wow  n  =  p(t)/p(t)  is  known  as  each  of  the  rays  x  =  x(t)  used  in 
the  construction  of  ^. -(t).   Consequently,  R.-  [^(t),  c(t)]  is  known 
for  all  t  >  t  at  each  point  of  jS.—{t).      The  problem  therefore  comes 
down  to  determining  how  e.-  (x)  varies  along  these  rays.   It  can  be 
shown   that  e.-  [^(t JJ  is  determined  from  the  following  equation: 


We  outline  a  proof  in  Appendix  I. 
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+       N 

2 


A.(x)U.(x)(e.(x)) 


lu.(x)(6.(x))^ 
1  •>   1  " 


X  =  x(t) 
i  =  s,  A,  f. 


r   t 


exp 


x=x(t  ) 
♦b  —   o 


T.  (t  )dt_ 


(3. IT) 


In  this  eq.uation,  x  =  x(t)  is  the  ray  through  x  =  x(t  )  associated 
with  the  fundamental  mode  under  discussion.   The  functions  T.(t),  i 


s^  A,  f  are 


T 


A 


2  (u  •  V  log  P) 


X  =  x(t) 


(3.18) 


T  = ^ 

i    2_^  ,  2  ,  2  ^2  -^ 
a  +  h  -4a  b 


P 


2  r 


^^  '    "^k     '  H  -  cj  [(n  •  V)u]  •  n   (3. 19) 


u, 
-^1 


2   2 
of-   tn  ^  2 


(cf  -  a^)Vlog  P  +^^-^-^V 


s,  f 


X  =  x(t), 


where 


+ 


a. 

1 


-  1;,    if  i  =  s, 
+  1,    if  i  =  f 


(3.20) 


+ 


A—  [^(t)j,  the  so-called  expansion  ratio  along  the  ray,  is  defined  as 

follows:   Let  dcr  denote  an  element  of  area  at  x  on  ><^(t  ).   Let  da.— 
o  .  -o    ^  ^  o'         1 

he  the  differential  area  on  «9. -(t)  which  corresponds  to  da  ,  the  cor- 

1  o 

respondence  being  effected  by  the  tube  of  rays  through  da  .   Then 

+ 

A.i  &(t)~]  ^p-~,  (3.21) 

^       -    da^ 

Since  U  =  p   (see  equation  (3.2)),  U  is  known  at  all  points  of  x  = 
X  (t). 
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+     +      +  2   + 
Noting  that  F.-  •  n.-  =  (e.-)  U.-  (see  equations  (3.1^)-)  and  (3.7) 

^  1    "^1      1     1 


(3.11))  and  combining  this  fact  with  equation  (3.I7)  we  find  that 

t 


+   +     + 
A.-  F.-  •  n.- 

1-1    -1 


X  =  x(t) 


»   +  -r,   +         + 

=  A.-  F.-  •  n- 

1  -^  1    ~i 


X  =  x(t  ) 


exp 


/T^(t)  dt 


t. 


(3.22) 


In  this  form,  we  see  that  the  normal  of  flux  of  wave  energy  out  of  a 


tube  of  rays  intercepting  the  area  dcr.-  on  the  phase  front  ^.-(t)  is 

proportional  to  the  influx  of  wave  energy  through  da  .   The  proportionality 
t  ° 

factor  exp  J_         T.  (t)  dt  ^  has  its  origin  in  coupling  of  the  fundamental 

to 
disturbance  waves  to  the  underlying  motion  of  the  medium  which  supplies 

energy  to  the  tube.   In  the  important  special  cases  where  u  =  0,  and 

where  both  P  and  u  are  constant  T. (t )  =  0  in  all  of  the  fuadamental 

■w  '1 

waves  and  equation  (3.17)  reduces  to 

\-   k(t)]  c.  [x(t),  n(t)]  (e.^  [x(t)])2  =  c(x^,nj(ei^^)2   ^^^^s) 

This  result  is  true  whether  H  is  constant  or  not.   If  in  addition,  H  is 

constant,  then  the  rays  are  straight  lines,  (see  Section  C  below), 

U  [x(t)]  =  p~  (t)  is  constant  along  x  =  x(t),  and  equation  (3.23)  reduces 

to 

\-  [^(t)]  (V  k(^2)^  =  (4,o^^  i  =  s,  A,  f  .       (3.2M 

In  this  case,  it  follows  from  the  theory  of  surfaces  that 


,iy,,2  .  !^C^-M 


(3.25) 


K.-  [x(t)] 

+ 
In  this  equation,  K^  [x(t  )]  and  K.~  [x("t )]  denote  the  Gaussian  curvatures 
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0       <^'  :^ 

at  x(t  )  and  x(t)  of  the  phase  fronts  ^^(t  )  and  /^ .    (t )  both  of  which 
are  assumed  to  have  curvature  at  these  points. 

Very  little  of  the  discussion  in  this  or  in  the  preceding 
sections  requires  modification  when  gravitational  effects  are  in- 
cluded.  Specifically,  the  addition  of  the  "gravitational"  terms  P  V  0 
and  SPV*  to  the  appropriate  members  of  equations  (2.1)  A  and  (2.2)5  A 
respectively  is  found  to  have  no  effect  on  the  formulas  for  the  disturbance 
speeds,  fimdamental  waves,  wave  energies  and  wave  fluxes.   The  presence 
of  the  gravitational  terms  is  felt  only  in  the  expressions  for  the  T^, 
i  =  s,  A,  f  given  above  and  then  only  in  the  addition  of  the  term 

V  <!>,      i  =  s,  f  .  (3.26) 


-  c  . 

T 


1 ,  gra V    h         2  2 
'   ^  c.  -  a  b 

1       n 


:  2   H     2  - 

b   -f  -  c.  n 
.  "  Hn    1  1 


to  T.  in  the  case  of  slow  and  fast  waves.   T.  requires  no  change  what- 
1  A 

soever. 


C.  Geometrical  Const2mction  of  the  Phase  Fronts 
Huygens'  Construction 
Again  let  b^ {t    )  be  a  small  neighborhood  of  the  point  x^  with 
normal  n  =  n(x  )  on  the  reference  phase  front.   Suppose  P,  H  and  u  are 
everywhere  constant.   Then  it  follows  immediately  from  equation  (3.7)  - 
(3.9)  that  on  each  ray  through  x  ^  E  =  P('t)/p(t)  is  a  constant  and  is 
equal  to  +  n  or  -n  according  as  the  ray  is  a  forward  or  backward  ray. 
Furthermore,  each  of  the  six  rays  through  x^  is  a  straight  line.   The 


""^The  fomiulas  mentioned  here  do  not  depend  explicitly  on  $.   However, 
the  parameters  r  and  b  appearing  in  these  formulas  are  of  course  affected 
by  the  presence  of  ^, 
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f± 
6^.  (t)^  i  =  s,  A,  f  are  thus  carried  parallel  to  themselves  along 

their  rays.   Explicit  formulas  can  be  given  for  the  phase  fronts 

evolving  out  of  ^  (t  )  as  soon  as  the  equation  for  jtoit    ),   x     = 

o  12 

f  (1^  |_)  is  specified  .   Here,  "we  shall  limit  ourselves  to  construct- 
ing these  phase  fronts  hy  p-urely  geometrical  means.   For  the  sake  of 
simplicity  it  will  he  assumed  that  u  vanishes  everywhere. 

In  Figure  3  we  show  the  position  of  6ji^   at  time  t  >  t  .   Wow 

according  to  equations  (3.3)  and  (3. 11) 
dx        dx     1 


n   •    —  =     n 


c[x(t),   n   (x(t)]  (3.27) 


dt  "O       dt         p(t) 

on  all  rays.      Since  H  and   P  are  constant  we  conclude  that 

dx 
n  ■    — -     =     c    (n   )     =     const.  ^  (3.28) 

dt  'O 

so  that,  in  particular,  for  the  fast  forward  ray 

n  •  —  =  c_  (n^)  =  const.  (3.29) 

'-   dt     f  vo 


The  vector 


— » 


P  F  =  c.  (n  )(t  -  t  )  n  (3.30) 

O         f   •'■O         0-0  V-J- J  / 


therefore  gives  the  displacement  of  5/^   along  the  direction  n  and 


the  plane  rj  (F)  drawn  perpendicular  to  P  F(P  )  at  F(P  ')  is  therefore 

tangent  to  Sv6'j,('t)  at  the  point  P  (x  =  x('t)).   This  fact  suggests  the 

following  construction.   At  each  point  P  of  g^{t    )   draw  the  vector 

P  F(P  )  =  c^(n  )(t  -  t  )n  .   Next  pass  a  plane"^  FfCP  )1  through  F(P  ) 
00      f"0        O'O  '-^i-^o-'  o 
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The  phase  fronts  evolving  out  of  an  initially  spherical  phase  front 


are  described  in  [lb] 


29  - 


which  is  normal  to  the  vector  P  F(P  ).   Then  clearly  the  envelope  of  these 


is  ^^(t). 


planes  is  ^^(t  j.   It  need  hardly  be  said  that  the  remaining  phase  fronts 
may  be  constructed  in  an  analogous  manner. 


Suppose  in  particular  that  ^ (t  )  is  a  "small"  sphere  of  radius 
e  located  at  the  origin  of  coordinates.   By  "small"  we  mean  that  e  is 
negligible  compared  with  the  distances  travelled  in  one  second  by  the 
^7(t)  evolving  out  of  ^  (t  ).   For  the  sake  of  definiteness  let  I£  be 
taken  along  the  x-axis  and  assume  that  the  ratio  r  -  a  /b  is  ,5.  The 
X-axis  is  evidently  an  axis  of  revolution  so  that  it  is  enough  to  con- 
sider the  cross  section  in  the  (x,y)-plane.   Such  a  cross  section  is 
shown  in  Figure  k.      The  ^J  "-curves  are  the  two  cusped  figures.   The 

h 

^.  -curves  are  two  points  (actually  two  spheres  of  radius  e)  located 
at  +  b;  in  this  type  of  wave  motion  5^  (t  )  splits  up  into  two  spheres 
congruent  to  oiJ  (t  )  which  propagate  along  the  H-direction,  with 
velocities  +  b.   The  29  "-curves  are  shown  as  one  oval-like  curve.  They 
are  actually  of  the  order  of  e  apart.   Figure  5  shows  why  this  is  the 
case.   Assuming  the  outward  normals  to^  (t )  correspond  to  the  forward 
propagating  waves,  we  see  thOj^  portion  of  the  wave  front  comes  from 
diametrically  opposite  points  on  the  sphere. 

The  curves  of  Figure  k   are  typical  of  those  associated  with  ■values 
of  r  in  the  range  r  <  1.   In  the  range  r  >  1,  the  Alfven  circles  are 
tangent  to  the  ^   curves  on  the  H-axis  but  the  general  appearance  of 
the  curves  remain  the  same.   The  case  r  =  1  is  somewhat  specialj  it  is 
discussed  in  detail  in  reference  [ibj . 
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Since  we  have  drawn  Figure  i|  as  if  6  were  zero,  we  may  thinly  of 
the  phase  fronts  pictured  there  as  emerging  from  a  point  source. 

These  phase  fronts  correspond  in  crystal  optics  to  I^bresnel's  ray  sur- 

12 
face.   We  refer  to  them  here  as  Friedrich' s  ray  surfaces   or  simply 

erf  t,  , 

as  spheres  or  wavelet s  and  henceforth  denote  hy  j.(t,P  ),    i   =  s,   A,t 

the  wavelets  (at  time  t)  that  have  issued  from  the  point  P  at  time 

t  <t/3 
o 

As  in  acoustics  and  geometrical  optics,  one  may  employ  these 
spheres  to  obtain  an  alternative  construction  of  the  phase  fronts 

r 

that  evolve  out  of  es^  (t  ).   This  construction,  known  as  Huygens ' 

construction,  consists  in  drawing  the  spheres  ^7  (t,P  )  at  each  point 

P  of  f^{t    )  and  then  forming  the  envelopes.   The  justification  of 

this  method  is  to  be  found  in  our  previous  discussion.   Consider,  for 

example,  the  construction  of  jaJ  (t )  which  was  explained  in  reference 

to  Figure  3.   Imagine  superimposed  on  this  diagram  the  sphere 

T"  ^(t,P  )  centered  at  P  on  Si^(t  )  where,  as  indicated,  the  normal 
'^f-'o  o        o 

is  ^  .   Consider  the  ray  on  ,^(t,P  )  whose  projection  on  n  is 
c_(n  )(t  -  t  )n.   Recalling  the  way  ,:!  (t,P  )  was  constructed,  one 
readily  verifies  that  this  ray  coincides  to  P  P  of  the  figure  and  that 
^j  (t,P  )  is  tangent  to  /$  „(t)  at  P.   From  these  remarks,  it  follows 
immediately  that  the  envelope  of  the  spheres  coincides  with  the  phase 
front  Ril^(t). 


We  saw  in  Sections  A  and  B  that  the  (unique)  ray  through  P  corre- 
sponding to  the  normal  n  to  5q^(t  )  at  P   (see  Figure  3)  is  the 


o 
con- 


12, 


They  were  first  described  by  K.  0.  Friedrichs  in  [l2] 
As  explained  above  j^j,(t,P  )  and  0*  „(t,P  )  coincide. 


31 


veyor  of  the  energy  of  the  associated  fundamental  wave.   In  view  of  the 

above  discussion  the  following  fact  is  apparent.   Only  those  rays  of  the 

spheres  jj".  (t,P  )  employed  in  Huygens '  construction  carry  energy  whose 

end  points  coincide  with  the  points  of  the  envelope  of  the  spheres  ^  . (t,P°), 

In  Figure  3^  P  P  is  one  such  energy-carrying  ray.   The  remaining  fast  rays 

issuing  from  P   (not  shown  in  the  figure)  carry  no  energy;  but,  of 

cotirse  they  help  in  the  visualization  of  the  envelope. 

The  constructions  that  have  been  described  in  this  section  are 
valid  even  if  H,  P  and  u  are  -variable,  provided  sufficiently  small  time 

intervals  are  chosen.   For,  the  basis  of  all  constructions  is  the 

dx 
relation  p  •  — *•  =  1  which  is  valid  whether  H,  P,  and  u  are  constant  or 
-   dt  ~         ■* 

not.   To  construct  a  phase  front  at  time  t  given  its  position  at  time 
t  ,  it  is  necessary  to  divide  the  time  interval  (t  ,t)  into  small  sub- 
intervals  and  starting  with  the  first  to  carry  out  the  foregoing  con- 
structions for  each  time  interval  in  succession. 

D.   Fermat ' s  Principle 
The  ray  equations  (3.7a)- (3. 9^)  all  have  the  form 

^  =  X  =  Vp'^Hx,?),    i  =  s.  A,  f  (3.31) 

dt 

where  Vp  is  the  gradient  in  p- space.   Let 

P  =  p(x,  x)  (3.32) 

be  the  solution  of  this  equation  for  the  components  of  p  in  terms  of 
the  components  of  x  and  x.   We  next  introduce  the  functions 

J'iCx,^)  = ?^7(x,  p  (x,  x) )  (3.33) 
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and  define  the  speeds  V.  as  the  solutions  of  the  equations 
+  r  + 


^i  K 1^ '  ^^  ^'    '  =  ^'  ^'  ^  ^^-^^^ 


+ 


V.  can  be  shovm  to  be  equal  to  the  magnitude  of  the  {+,    i)-th  ray 

velocity.   More  important  here  is  the  fact  that  V.  =  V(x,x)  is,  by 

definition,  a  homogeneous  function  of  the  degree  zero  in  the  components 

of  X.   It  follows  that  a  change  of  parameter  leaves  the  functional 

+ 
form  of  V.  unchanged  if  by  x  we  always  understand  the  derivative  with 

respect  to  the  param.eter. 

An  {+,    i)-ray  may  now  be  characterized  as  the  curve  which  makes 

the  time  of  travel 

VJo-  il'  =  J        ^,  (3.36) 

an  extremum.   All  competing  curves  in  this  variational  principle  are 
required  to  pass  through  the  points  x  and  x.,  when  s  =  s  and  s., 
respectively;  the  curve  parameter  s  is  the  arc  len^h.   This  is  Fermat's 
principle.   We  refer  the  reader  to  R.  K.  Luneberg's  notes  [ii]  for  the 
proof  of  this  principle  and  a  discussion  of  its  relation  to  other 
m.ethods  of  characterizing  the  rays.   Here,  we  merely  note  that  the 
application  of  Fermat's  principle  is  limited  in  practice  to  those  cases 
where  it  is  easy  to  express  p  in  terms  of  x  and  x  and  to  find  the  roots 
of  equation  (3.35). 
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rV.   PROPAGATION  OF  ALFVEN  AND   FAST  FRONTS  IN  INHOMOGENEOUS  MEDIA 


A.   General  Remarks 


To  illustrate  the  theory  of  Sections  II  and  III,  we  consider 

the  propagation  of  Alfven  phase  fronts  in  a  stationary  but  in- 

+ 
homogeneous  medium.   In  this  case,  the  ray  and  e.- strength  equations 

can  he  written  in  the  form  [see  equations  (S.T)  and  (3.23)] 


dT     -  N 


V 


A  ' 


(Kl) 


to     t 
^  A    A  A 


x(t) 


=  ^^a')'=a 


x(t  ) 
•»  o 


(^.2) 


Here,  N  is  any  conveniently  chosen  normalization  constant,  and  the 
new  parameter  T  is  defined  by  the  equation 

T 


t-t^   =   ^   /   /P[X(T')]  ^-'. 


{k.3) 


±  ±  ± 

A  ,  it  will  be  recalled, is  the  ratio  dcr./da  of  the  areas  da.  (t)  and 


A' 


A'   o 


A' 


da  =  da(T  )  on  the  phase  fronts  4/v(t)  and  ^(t  ),  respectively,  cut 

out  by  a  tube  of  rays  through  da  , 

°  +  -  ■*"  - 

According  to  equation  (4.2),  e    (s  )  depends  both  on  k  {k   )   and 

on  the  direction  of  the  wavefront  normal  n(since  e.  =  (WTP  )^).  Thus, 

■^  An  ' 

it  appears  that  the  determination  of  e     requires  not   only  the   solution 

of   {k.±)  but  also  of  the  system  (3.9'b).      However,    it  can  be  shown  that 

+ 


d  a. 


d  a. 


G       P' 
A 


x(t) 


c  P2 

^A 


X  =  x(t^) 


(ii.M 


3h  - 


so  that   eqiiation   (4.2)   reduces  to  the   surprisingly  simple  relation 

+  -      1 


+  o     i 


=      (e-)2p2 


=  x(t) 


(h.-j) 


X   =  x(t    ) 


Thus  the  e~  actually  depends  only  on  the  local  density  and  are  in- 
dependent of  the  expansion  ratio.   Moreover,  if  P  is  everywhere 
constant,  €  and  g"  remain  constant  along  their  rays. 

To  prove  equation  (4.5),  we  integrate  V  •  v  over  the  volume 

»• 

AV     of  the  tube  of  rays  between  dcr      (t)  and  da   (t   +  A  t).      Next, 

+ 
we  employ  the  divergence  theorem,    divide  both   sides  by  A  V     and 

-±  ±       /  -  -     r 

make  use  of  the  relations  v„      =  v.    dt/dT, v  .    .    n  =  U.        see  equations 

(4.1),    (3.11)  and   (4.3)3.      Letting  At  approach  zero,   we  find  that 

(cf.    [4]    p.    87). 

+ 


V  •  V 


=     ^—  log 
dT        ^ 


+  i 

da^  %      /P^2 


L^%%o\g/  J 


(4.6) 


Here   P    and  U.   „  denote  the  density  and  normal  speed  on  aij  (x    ).      It 
■^  o  A,  0  o 

may  be  remarked  that  this  resiilt   is  completely  general;   u  need  not 

vanish.      In  addition,    it  applieSjWith     obvious  change     of  notation, 

to  slow  and  fast  waves.      Equation   (4.4)  now  follows   from  the   fact 

that  V  •   V?  =   +  E'^  V  •    H  =  0. 
♦"A        —  *«' 


B.   Propagation  of  Alfven  Fronts  in  a  Dipole  Field 

We  turn  now  to  the  problem  of  determining  the  propagation  of 
Alfven  fronts  in  a  dipole  field.   In  this  case,  the  basic  equations 
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are: 


H     = 

-   V 

(M  cos  0/R^   ), 

dR 
dT 

-    2 

M  cos  0 
N     r3 

Rd  e 

M  sin  e 

d  T       N  r3 


(h.8) 
ih.9) 


i^o)  = 


^o  =     ^  ^^o^  *o^  ^^-1°) 


Here,  M  is  the  magnetic  moment  of  the  dipole,  R  is  the  distance 

from  the  origin  and  0  is  the  polar  angle  measured  from  the  vertical 

dipole  axis.   Equation  (U.IO)  is  the  equation  of  the  reference  surface 

in  polar  coordinates  (t)  being  the  azmuthal  angle.   For  the  sake  of 

'o 

simplicity,  we  consider  only  the  rays  anti-parallel to  the  field  -  i.e. 

those  which  are  directed  into  the  north  magnetic  pole. 

Eliminating  t  between  equations  (^.8)  and  {h.9),   we  find  the 

equation, 

R  =  k  sin^  e,  (U.n) 

o 

,^  __  ^    __  li!o^  (,.,„ 

sin  d  sin  9 

o  o 


Ji 


for  the  ray  through  the  point  (R  ,  9  ,  ({)  )  on  j^(t  ). 

We  now  use  equation  (4.11)  to  eliminate  9   in  equation  (4.8). 
On  integrating,  we  obtain  the  result 

R  h         R/k_, 

T  -  T   =  —  /   (+)  ■  =  /    (+)  f=   dx  .        a   13) 

°   2M  {     x/l  -  (r/kj   2  M  ^^  .    X/II^  ^   ^^ 

o  o'  o 


36  - 


The  minus  sign  in  the  last  member  is  to  he   applied  over  that  range 
of  R  in  which  0  <  0  <  jt/2  and  the  plus  sign  in  the  0-range 
n/2  <  9  <  -a.    Setting 

o 
where  (see  reference  [2^]  ) 

B  =  /    ""    "^"^        =     .91^28  (4.15) 

and 

W     =       2M/B,  (4.16) 

we  find  that   equation   (4.13)   reduces  to 

T-T       =    k^  [I   (R/k^)  -  KR^kJ]      ,  (4.17) 


where,  in  virtue  of  equation  (4.3)^ 

^  -  ^o  =  3=   J     Vp[(t')]  ^-'  .  (^-18) 


^o 

l(x)  will  he  recognized  as  an  incomplete  gamma  function.   Since 
these  functions  have  been  tabulated  (see,  for  example,  reference  [24]  ) 
equations  (4.17)  and  (4.18)  furnish  the  solution  to  our  problem. 


For  simplicity,  suppose  now  that  ^{r    )   is  the  plane  whose 


equation  is 

R    =  d/(sin  9     cos  (t)  )  (4.19) 

o       '       O       'O        1 
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d  "being  the  distance  fi-om  the  dipole  in  the  equatorial  plane.   If, 
in  addition,  we  assume  that  P  is  everywhere  constant  then 

t  -  t   .  2M   ^^  _    ^^  ^j^^20) 

o     2         o 

and  T  is  proportional  to  t.   In  Figure  6,   we  have  pictured  the 
evolution  of  the  phase  fronts  in  the  meridianal  plane  ^     =0  from 


..........     .  .  ^ ..        _       j{(Q^      .^ 


the  time  t  =0  onward.   The  normal  n  to  ji4(0)  is  assumed  to  point 
o  — o    •^  -^ 


to  the  left  and  to  define  the  forward  direction.   The  curves  shown 
correspond  to  the  t- values  0,  .1,  .2,    .3  and  .k.      To  keep  the  diagram 


imcluttered,  we  show  ,<j.  (t)  on  the  upper  half  of  the  diagram  and 
i^^{'^)   on  the  lower  half.   The  complete  picture  is  obtained  by  adding 
to  the  curves  already  present  their  reflections  in  the  (x/d)-axis. 
Since  the  rays  associated  with  the  Alfven  mode  are  the  field  lines  and 
since  these  lines  converge  at  the  origin,  one  might  expect  the  phase 
front  to  develop  a  sharp  peak  as  it  nears  the  dipole.   We  see  from 
the  figure  that  this  is  indeed  the  case.   The  dashed  line  is  the  field 
line  on  which  the  peak  develops.   Its  equation,  obtained  from  the 
figure,  is 

R  =  1.15^1  sin^  Q.  {h.2l) 

It  may  be  imagined  that  the  formation  of  the  peak  corresponds  to 
a  concentration  of  energy  at  the  peak.   However,  this  is  not  the  case, 
since  according  to  equation  (^.5)  "the  disturbance  strength  remains 
constant  along  a  ray. 

We  have  already  established  the  close  relationship  between  the 
time  harmonic  theory  and  the  pulse  theory.   This  relationship  allows 
as  to  regard  to  the  ^(O)  in  Figure  6  either  as  an  initial  surface 
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of  discontinuity  or  as  a  fixed  sirrface  in  space  carrying  a  disturbance 
vhose  strength  varies  harmonically  in  time.   In  either  case,  if  Figure 
6  is  to  make  sense,  it  must  be  assumed  that  the  disturbance  at  ^^  (O) 
is  of  the  form  eR     +^fln^An  ^^®®  equations  (2.28))  and 
that  £ ft  c  f  0.   If  for  example  e.  ^  then  the  backward  propagating 
front  ^^^(t)  must  be  ignored  since  it  carries  no  disturbance. 

An  "earth"  of  radius  ad,  0  <  a  <  1,  can  be  indicated  on  Figure 
6  by  drawing  a  circle  of  radius  a  d  centered  at  the  origin.   The 
semicircle  C  in  Figure  6  corresponds  to  an  earth  of  radius  l/5  d. 
Let  us  for  the  moment  regard  ^(O)  as  a  surface  of  discontinuity. 
Then  J1a(t)  and  jj<>.  (f)  are  Alfven  discontinuity  fronts.   The  figure 
indicates  that  jQ    (t)   reaches  the  earth  first  at  the  point  P  where 
the  peak  just  makes  contact  with  C.   Then  the  disturbance  spreads 
over  higher  and  lower  latitudes  as  the  front  behind  the  peak  passes 
over  C.   (We  are  ignoring  here  the  effect  of  the  earth  on  the  dis- 
turbance wave).   If  we  imagine  that  the  region  to  the  left  of  ^(O) 
is  undisturbed  while  the  region  to  the  right  is  disturbed  then  we 
may  safely  conjecture  that  the  disturbance  at  each  latitude  will 
undergo  two  rapid  changes.   At  P,  for  example,  the  first  change  occurs 
when  the  peak  reaches  the  earth  and  the  second  when  the  disturbance  on 
the  surface  element  located  at  P'  has  reached  P  after  traversing  the 
dashed  field  line  in  the  counterclockwise  direction. 
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C.   Propagation  of  Pkst  Waves  in  the  Earth's  Magnetosphere 

In  the  earth's  magneto  sphere,  the  soimd.  speed  is  much  less  than 

the  AlfVen  speed.   For  this  reason,  we  may  ignore  all  but  the  first 

2   ? 
terms  in  expansions  involArLng  the  parameter  r  =  a  /b  .   We  then  find 

that  the  slow  disturbance  speed  vanishes  and  that  the  fast  disturbance 

speed  and  fast  ray  equations  reduce  to 

.2  -1 


c^  =  (n  E""   p"^)  =b,  (4.22) 


il  =  b^p  =b  n,  (lf.23) 

dt        •-     * 

^  =  -  p  Vb  .  (k.2k) 

d  t 

-1 
Here,  we  have  used  the  fact  that  p  =  c  . 

For  simplicity,  we  henceforth  limit  the  analysis  to  magnetic 
fields  whose  magnitude  may  be  expressed  in  the  form[Hl=  const,  R   , 
where  R  =x/„2  ^  2.   This  choice  allows  us  to  represent  a  uniform 
field  (n  =  O),  the  field  of  a  line  current  (n  =  l),  a  two-dimensional 
dipole  (n  =  2)  and  the  field  in  the  equatorial  plane  of  a  three  di- 
mensional dipole  (n  =  3).   Note  that  any  two-dimensional  multipole 

field  has  the  desired  form.   This  follows  the  fact  that  the  complex 

*  *   „   -n 

conjugate  H  =  H  -iH  of  such  a  field  has  the  form  H  =  A  z  , 
u-to         X   y 

z  =  X  +  iy,  from  which  it  follows  that  IT  =  H  H  =  |A|   R 

Eliminating  p  between  equations  (il-.23)  and  (if. 2i4-),  we  find  that 

J:(\    tl]    =  2vb-2  .  (1^.25) 

dt  \  b'^  d  t/ 


-  k-0 


Introducing  the  polar  coordinates  (R^,  ^)   and  assuming  that  the  density 
depends  solely  on  R,  we  conclude  from  this  equation  that 

2  • 

^  i  =  ^^  =  MR,,  ^,;nj,        ik.26) 

o 

where  the  initial  value  (t)   is  to  calculated  from  equation  {k.23). 

'  o 

Next^  we  dot  equation  (4.23)  with  itself  and  and  find  that 

R   +  R^  <{)  =  t^.  (4.27) 

The  last  two  equations  lead  directly  to  the  differential  equations 

^2  „2n  +  1 
dt       PR  _     ^  (1^^28) 


^«   "  CC  /r2-  -^  2_  ^2^^  ^.1 


a.  Lop2 


Li.  ~-      ""°^  (4.29) 

<iR     \/^^^.a^L%-^      ' 


Here,  we  have  made  use  of  the  fact  that 

J, 
h  =  a  p"^r""  ,  ■  (4.30) 

a  "being  same  positive  constant.   When  p  is  everywhere  constant  we 
obtain  the  explicit  results 


(4.31) 


t 

-  t 
0 

(n+l)a 

J  L  2n+  2         2     2-1 
WR               -Lap 
0 

-     /„2n+  2     -,2     2     -1 
/^o          -   ^o  «     P 
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r"""^  COS  [(n  +1)  (())  +  P)J   =  «V\/^  ^^-32) 

In  equation  (4.32)  p  is  a  constant  of  integration  which  is  determined 
by  evaluating  the  lefl:  member  of  (R  ,<()  )  on  the  phase  front. 

In  the  magnetosphere  p  is  a  complicated  function  of  R  and  numerical 
Integration  is  necessary.   Employing  the  density  profile  in  the  magneto- 
shere  W.  E.  Frances  et  al,  [^25]  have  calculated  several  ray  paths  in 
the  equatorial  plane  issuing  from  a  point  source.   These  authors  obtain 
their  ray  equations  by  postulating  Femiat's  principle  using  b  (see 
equations  (4.22)  and  (3.36))  as  the  ray  speed.   It  would  be  worthwhile 
to  extend  these  calculations,  now  that  they  have  been  placed  in  their 
proper  context,  to  derive  the  strengths  and  form  of  the  fronts  that 
originate  in  suitable  extended  sources.   It  is  our  intention  to  in- 
vestigate this  question  in  a  future  report. 
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V.   THE  LAWS  OF  REH^CTIOW  AND  REIEACTION 

A.   Preliminary  Remarks.   Snell's  Laws 

It  has  thus  far  been  assiomed  that  the  undisturhed  medium  is  every- 
where continuous.   Suppose  now  that  a  surface  of  discontinuity  ciC7 
separates  two  different  states  of  the  undisturbed  medium.   For  the  sake 
of  simplicity,  we  shall  require^C/  to  be  fixed  in  space  although  our 
results  apply,  with  obvious  modifications,  to  the  case  whereof  is  in 
motion  and  changes  shape.   Imagine  that  one  of  the  three  fundamental 
disturbance  waves  is  incident  on  oO  .   Denote  by  UJ  . ,   the  incident 
phase  front,  where  i  =  s.  A,  f,  and  suppose  that  the  position  of  w  . 
and  the  disturbance  strength  at  each  point  of  WV  .  are  known  at  the 
reference  time  t  .   The  basic  problems  of  the  theory  are  then  to  deter- 
mine (l)  the  laws  of  reflection  and  refraction  and  the  motion  of  the 
scattered  phase  fronts,  (2)  the  strengths  of  the  scattered  waves  in 
terms  of  the  strength  of  the  incident  wave. 

We  shall  be  concerned  here  with  the  first  of  these  problems. 
Formulas  for  the  strength  of  the  scattered  waves  emerging  from  various 
discontinuities  are  given  in  references  [l]  [22]  and  [23^  . 

Unless  otherwise  mentioned,  it  is  assumed  that  the  state  of  the 

unperturbed  medium  on  each  side  of  the  surface  of  discontinuity  o^ 

is  constant  and  that  the  incident  phase  fronts  are  planar.   It  should 

therefore  be  stressed  that  all  results  obtained  for  this  special  case 

apply  locally  (in  the  neighborhood  of  the  curve  of  intersection  of  c*-? 

with  the  incident  front  2C/.  )  when  tJ-^  or  the  incident  wave  is  curved 

1 

or  when  the  luidisturbed  media  separated  byc>^are  non- uniform. 
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Figures  7  and  8  depict  the  basic  geometrical  set  up.  "^represents 
the  sui'face  of  discontinuity  separating  the  upper  half-space  6^'  from 
the  lower  half-space/\;|_.   The  vector  H  is  the  prevailing  magnetic  field 
at  the  origin  0.   The  local  plane  of  incidence  at  the  origin  0,  deter- 
mined by  the  normal  on  the  incoming  wave  front  and  the  unit  normal  N 
on  ct'  at  0,  is  assumed  to  coincide  with  the  (x,z)-plane  and  is  denoted 
by  I-J-,  As  shown,  N  is  directed  along  the  positive  z-axis.   The  vector 
H  is  the  projection  of  H  on  the  plane  of  incidence.   The  quantities 

Tp^,    0  <  f^  <  ■n/2   and  0„,    -rt  <  0„  <  tc,  are  the  angles  between  H  and 
ti  —   n  —  il         xl  <*• 

H  and  H  and  N,  respectively.   We  shall  have  occasion  to  refer  to 
the  angles  P  and  7  of  Figure  7  later. 

In  the  above  remarks  and  in  the  drawing  of  Figure  J,    it  has  been 
implicitly  assumed  that  H  is  continuous  across c>c7.   This  need  not  be 
the  case;  ^J-'ma.y,    for  example,  be  a  stationary  hydromagnetic  shock. 
Accordingly,  we  shall  write  H  for  the  magnetic  intensity  in  t<  and  use 
H'  to  denote  this  quantity  in  6l_' .      In  general,  if  Q  is  any  quantity 
in  ^  we  employ  Q'  to  denote  the  corresponding  quantity  in  ^^'  whenever 

Q  4  Q'- 

Figure  8(a)  shows  the  phase  fronts  in  the  plane  of  incidence.  ^i(/ . 

denotes  the  incoming  phase  front,  propagating  toward  <0  from  below.   We 

take  l^.  to  be   fast,  Alfven  or  slow  according  as  i  =  1,  2,  or  3.  w±2_> 

U/. ^  and  u/1  denote,  in  the  order  stated,  the  reflected  fast,  Alfven 
i2       13  

and  slow  phase  fronts  and  iP-i^,    LV-c-   and  UU ■  r   the  refracted  slow,  Alfven 
and  fast  phase  fronts  -  all  produced  by  the  incident  wave  of  type  (i). 
Figure  8(b)  shows  the  disposition  of  the  normals  n.  and  n.  .  on  these 
wave  fronts;  to  keep  the  diagram  uncluttered  only  three  normals  have 
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actually  teen  drawn.   The  normals  point  along  the  direction  of 

propagation  of  their  fronts  (see  equations  (2.9)  and  (2.11)).   The 

angle  9.,  -n/2  <  9.  <  rt/2,  is  the  angle  between  N  and  the  incident 

phase  normal  n..   Here,  and  in  the  following,  angles  in  the  plane 

of  incidence  are  regarded  positive  or  negative  according  as  their 

sense  is  clockwise  or  counterclockwise.   We  accept  at  this  point 

the  easily  established  fact  [22]  that  the  phase  normals  n.  .  of  the 

scattered  waves  all  lie  in  the  plane  of  incidence.   Then  we  denote 

^y  ^-  -^  J  =1;  '•')    6,  angles  (in  the  plane  of  incidence)  between 

N  and  the  j-th  scattered  wave.   Specifically,  when  j  =  1,2  or  3, 

the  angle  0.  .,  -  -  <  9.  .  <  -,    is  the  angle  between  N  and  the  j-th 
ij    2  -  ij  -  2' 

reflected  wave.   When  j  =  ^,  5,  or  6,  6.  .     is  the  angle  between  N  and 
the  j-th  transmitted  wave. 

In  the  plane  at  incidence  I  (see  Figure  7)  the  incident  and 
scattered  phase  fronts  intersect  in  a  common  point  0  which,  it  will 
be  supposed,  moves  with  the  velocity  V  =  V  x  at  the  instant  depicted. 
V  is  called  the  trace  velocity.  Let  U.  >  0  be  the  speed  at  0  in  the 

,<v  1 

normal  direction  n.  (see  equation  (2.10))  of  the  incident  phase  front 

and  let  U.  .  denote  the  corresponding  normal  speeds  of  the  scattered 

ij 

waves.   From  the  fact  that  *u/.  and  the  ^//T  .  intersect  at  0  and  move 

together  with  common  velocity  Y,   we  conclude  immediately  that  for  any 

i,  i  =  1,2,3,  that 

V  •  n.   =  U.,  (5.1) 

—   ~i      1 


and 


V  •  n.  .  =  U.  .     j  =  1,2,. ..,6.         (5.2) 
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We  have  seen  that  normal  speeds  are  direct ion- dependent.   This 
dependence  on  direction  is  made  explicit  in  the  following  manner: 


where 


U.   =  u-r  •  n.  +  c. 

1     ^I   -1  -   1 


b 


i  =  1.  2,  3, 


1  .  2 


1+  r  _^  1 
2     2 


(l+  r)  -  kv   cos  t{/^   cos  {9  -9    ) 


-.3-  ^ 


H 


1  H' 


(5.3) 


b.h) 


=     cos  Tprr     I  COS   ( 


H 


)  I  . 


(5.5) 


b 


=  < 


1  +  r    1 
-~2      2 


1  2 


(1+  r)2  -  ^r  cos2  fjj_   cos^CS  -fl  ) 

n  - 


(5.6) 


In  equation  (5.3)^  u  is  the  projection  of  u  on  the  plane  of  incidence. 
The  quantities  c  ,    Cp  and  c  are,  in  the  order  stated,  the  slow,  Alfven 
and  fast  disturbance  speeds.   The  quantity  r,    it  will  be  recalled,  is 
the  squared  ratio  a  /b  of  the  sound  speed  a  to  the  Alfven  disturbance 


speed  b  =  \/|a  H^p-L   Equation  (5.3)  follows  directly  from  the  definition 
of  U.  (see  equation  (2.17 )).   Similarly,  the  expressions  for  the  c. 
follow  directly  from  the  definition  for  the  Alfven,  slow  and  fast  dis- 
turbance speeds  (see  equations  (2.I9)- (2.2l) )  on  substituting  for  n 

J 

and  H  the  expressions 

cos  0.  N  +  sin  e.  X   ,  (5.7) 

1    ^  1   *o 


n . 


3    -  J  *o' 

H     =     H  cos  i}J^   (cos  9^  N  +  sin  9     x)    +  H  sin  ?//„  z.  (5.8) 
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The  corresponding  formulas  for  the  U.  .  and  c.  .  may  be  obtained  from 
equations  (5.^)- (5.6)  by  the  simple  expedient  of  replacing  v,   b,  ijj^ 
and  0  by  r',  b',  ^  ,  and  Q„x   and  replacing  6.   by  the  appropriate  Q.  ., 
j  =  1,2,  ...  6.   In  accordance  with  oiir  previous  notation,  the  indices 
j  =  1,2  and  3  si's  to  be  associated  with  the  reflected  fast,  Alfven  and 
slow  fronts  respectively;  and  the  indices  j  =  ^,  5^  and  6  with  the 
slow,  Alfven  and  fasx  transmitted  waves. 

Combining  equation  (4.1)  and  (4.2)  with  (4.3)  and  the  corresponding 
formula  for  the  scattered  waves  we  find  for  each  i,  i  =  1,2,3,  that 

(V-  uj)  •  n.   =  ic.,  (5.9) 

(V-u).n.  .   =+c   ,    j  =  1,2,3,  (5.10) 

(V  -  u-  )  •  n   =  +  c      j  =  3,^,5.  (5.11) 


Sett  ing 


and 


V  =   I  V  -  Uj  I,  (5.12) 


V   =  I  V-  u'J,  (5.13) 

and  denoting  by  o;  and  a'  the  angles  that  the  velocities  V  —  u-p  and 

V  -  u'-r  make  with  the  surface  normal  at  N,  we  arrive  at  the  hydroraagnetic 

analogue  of  Snell's  laws;  namely 

IcosCe.  .  -  a)|     I  cos  {6.-  a)    \     i 

^-^ = =  —  ,  j  =  1.2,3,    (5.14) 

V 


—  '  J  =  ^>5,6.    (5.15) 
V 


^id 

c. 

X 

cos(0ij   ■ 

-  c^ 

cos(0i-  a') 

'^ij 

=1 
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Introducing  into  these  equations,   the  quantities 

P.     =    ^  .  (5.16) 

1  Cj_/b 


c^j/t,  J   =  1,2,3 

i.AS  j  =  3A,5. 


^ij=        '^,  (5.17) 


c 


and  the  Alfven  numbers 


M,    =  -~-  ,  (5.18) 

^     (v/b) 

Mb'  =  -3T-  ^  (5.19) 

^      V'/b' 

we  may  write  Snell's  laws  in  the  following  useful  form: 

p. .  I  cos  (0. .  -  a)  I  =  p.  Icos  (e.  -  a)  I  =  — ,        (5.20) 

p.  .  I  cos  (9.  .  -  a')  I  =  p.  Icos  {e.-a  '  )|  =  —  .       (5.21) 

b 

Regarding  (p.  .,Q.   .)  as  the  polar  coordinates  of  points  in  the  (p,6)- 

plane,  p  being  the  radial  variable,  we  see  immediately  that  the 

(p. .,  9.   .),  j  =  1,2,3,  all  be  on  a  line  through  (p. ,9.)  whose  equation 


is 


Z:      p  I  cos  (0  -  a)  I  =  i-  =  p^  |cos  (9.  -  cc)    \  (5.22) 

b 

while  the  (p.  .,  0.   .),  j  =  4,5,6  all  be  on  a  line  i'  through  (p. ,9.) 
whose  equation  is 

i':  p  I  cos  (9  -  a')  =  =  p.  Icos  (9.  -  a' ) I  .  (5.23) 
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Suppose  now  we  plot  p.  .j,  j  =  1,2,3  against  6.  .  on  the  same  (p,0) 
polar.  The  result  is  a  family  of  five  curves  -  the  Alfven  and  slow 
curves  consist  of  two  branches  each.   It  follows  immediately  from 
equations  (5.20)  and  (5.22)  that  the  points  of  intersection  of  the 
line  i   with  the  curves  of  this  family  are  solutions  of  equation  (5.20) 
and  hence  fiu-nish  the  angles  of  reflection  and  the  reciprocal  disturbance 
speeds  of  the  reflected  waves.   These  remarks  are  illustrated  in  Fig- 
ure 9(a).   The  parameters  x,    Q„,  i\J^   and  a  are  .5,  ^5^0  and  0 
respectively.   Since  a  =  0  the  line  I   is  vertical.   The  distance  of 
I   from  the  vertical  W-axis  is  evidently  mT  .   Figure  9(^)  depicts  the 
corresponding  diagram  for  the  transmitted  waves. 

We  shall  discuss  our  graphical  method  at  greater  length  later. 
First,  however,  we  wish  to  examine  its  relationship  to  Huygens' 
classical  construction. 

B.   Huygens'  Construction  of  the  Scattered  Phase  Fronts 

We  assume  that  the  mediiwi  in  U^  is  at  rest  (u  =  O)  and  consider 
the  construction  of  the  reflected  fronts  by  Huygens'  method  which, 
in  brief,  is  the  following:  At  the  instant  the  incident  phase  front 
reaches  the  point  0  of  the  interface  (see  Figure  7,  8),  0  becomes  an 
induced  point  source  and  emits  a  wavelet  or  sphere  (see  Section  II  C). 
To  construct  the  reflected  phase  fronts  at  time  t'  one  simply  con- 
structs the  envelopes  of  the  spheres  emitted  at  all  times  t  <  t'. 

In  Figure  10,  0'  represents  the  trace  point  at  time  t'  and  0  is 
the  position  b"   seconds  earlier.   It  is  assumed  that  0'  moves  with 
the  constant  velocity  V  =  V  x  and  that  V  >  0.   The  distance  of  0' 
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from  0  is  therefore  V/b.   Taking  r,  ip     and  0„  to  be  2  ,  0  and  ^5 
respectively,  the  sphere  through  0  is  the  sirrface  of  revolution  ob- 
tained by  revolving  the  f , ,   A  and  s  -curves  of  the  I^lgure  10  about 
the  H-axis.   The  UJ.  in  Figure  10  represents  the  three  possible  types 

of  incident  phase  fronts  and  the  lAj.  .   represent  the  associated  family 

-^  J 

of  reflected  phase  fronts.   The  "(/.  .  and  the  dashed  extensions  of  the 
{jf.    are  tangent  to  the  f  ,   A  and  s  -curves,  since  by  construction 
these  phase  fronts  must  be  the  envelopes  of  the  wavelets  emitted 
earlier  and  up  to  b   seconds  later  than  the  sphere  shown.   It  is 
perhaps  unnecessary  to  stress  that  the  hemisphere  in  the  transmission 
region  ^'  indicates  merely  where  the  disturbance  would  have  been 
were  there  no  discontinuity  across  ^tC/.   The  hemisphere  through  0 
appropriate  to  ^',  will  in  general,  differ  from  that  shown  and  as 
in  geometrical  optics  may  be  employed  in  the  construction  of  the 
transmitted  phase  fronts. 

The  geometry  of  Figure  10  changes  only  in  an  obvious  manner  when 
the  parameters  r,    0„  and  V/b  are  varied.   Thus  the  Alfven  spheres  are 
tangent  to  the  s  or  f^  curves  according  as  r  >  1  or  r  <  1;  but  other- 
wise varying  r  effects  only  the  detailed  features  of  the  figure.  Vary- 
ing 9     has  the  effect  of  rotating  the  f^ ,  A  and  s_  curves  rigidly 

n  1    d       j 

about  0,   Finally,  varying  V/b  moves  0'  closer  to  or  farther  from  0 

15 
and  governs  what  may  be  termed  the  angles  of  "critical  reflection"   . 


15 

A  reflected  front  of  type  j  due  to  a  incident  wave  of  type  i  is  said 

to  be  critically  reflected  when  0.  .   just  becomes  equal  to  +  9O  .  The 

angle  9.    at  which  this  occurs  is  called  the  angle  of  critical  reflection 

for  the  j-wave. 
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For  example,  if  0'  falls  with  the  f  -oval,  a  fast  phase  front  is  not 
reflected  since  there  are  no  tangents  to  the  f, -curve.   In  addition, 
ve  iray  mention  an  anomalous  sort  of  reflection  which  occurs  when  0' 
coincides  with  0"  of  Figure  10.   In  this  case,  it  is  easy  to  see 
that  an  incident  slow  front  produces  reflected  slow  and  Alfven  fronts 
which  lie  in  the  same  quadrant  as  the  incident  wave  itself.   The  flux 
of  energy  (see  equation  (3.1^))  associated  with  each  of  these  fronts 
is,  however,  in  the  proper  direction,  that  is,  away  from  the  inter- 
face 

It  has  thus  far  been  assumed  that  H  lies  in  the  plane  of  in- 
cidence  i.e.,  that  Tp„  =   0.   Increasing  ip^   away  from  zero  has  the 
effect  of  turning  the  H-axis,  the  axis  of  revolution  of  the  wavelets, 
out  of  the  plane  of  the  page.   Tlie  points  of  tangency  with  the  wavelets 
of  the  incident  and  scattered  fronts  also  move  out  of  the  plane  of  the 
page  and  are  no  longer  coplanar.   Although  Huygens '  method  seems  now 
to  require  three  dimensional  figures,  it  is  nevertheless  still  possible 
to  perform  the  construction  by  the  procedure  described  above.   One  need 
only  work  with  the  projection  onto  the  plane  of  incidence  of  the  locus 
of  points  on  the  wavelets  at  which  the  normal  n.  is  parallel  to  the 
plane  of  incidence.   The  parametric  equations  of  this  locus  may  be  ob- 
tained from  the  ray  equations  (3.T)-(3.9)  "by  substituting  the  expressions 
for  n  and  H  given  in  equations  (5.7)  and  (5.8).   Figure  11  depicts  this 
locus  for  various  values  of  ip  ;   r,  and  9     are  again  assumed  to  be  l/2 
and  ^5  .   The  Alfven  locii  are  not  shown. 
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To  explain  the  essential  elements  of  the  relationship  of  Huygens ' 
construction  to  the  graphical  procedure  of  Section  A,  it  will  suffice 
to  show  how  to  obtain  Figure  9  a  from  figure  10;  in  both  figures 
r  =  1/2^  d     =  k^     and  Tp„   =  0.   We  maintain  that  figure  9a  may  be  ob- 
tained from  figure  10  by  l)  constructing  the  pedal  curves  to  the 
curves  of  figure  10,  and  then  inverting  these  pedal  cmrves   with 
respect  to  the  unit  circle  centered  at  0. 

To  justify  this  statement  we  note  that  the  pedals  2.    from  0  to 
the  curves  f  ,  A  ,    s     are  2.    =  b~  c.(n. )n.,  i  =  1,2,3,  n.  being  the 
normal  to  the  i-th  curve.   This  fact  follows  from  the  way  the  wave- 
lets are  constructed  and  the  fact  that  figure  10  corresponds  to  an 

elapsed  time  b    (=  t'-t)   seconds.   The  inverse  point  of  2.   with 

1 

respect  to  the  unit  circle  at  0  is  evidently  the  vector  directed 
along  n  of  length 

P.  =  —} — rr  '   i  =  1^2,3  (5.21+) 

1     Ci(ni)/b 

Comparing  this  equation  with  equation  (5.I6),  we  conclude  with  aid 
of  equations  (5.7)  and  (5.8)  that  the -corresponding  (p,0)-polars 
and  are  in  fact  those  shown  in  figure  9{b.). 


The  pedal  curve  of  a  curve  C,  whose  equation  is  x  =  x(|),  is  con- 
structed as  follows:  Let  T(4)  be  the  tangent  to  C.   From  the  origin 
0,  draw  the  perpendicular  to  T(|)  and  let  y  =  y(|)  denote  the  vector 
from  0  to  the  foot  of  this  perpendicular.  The  c-urve  traced  out  by  y 
as  4  varies  over  its  domain  of  definition  is  the  pedal  curve  of  ^. 
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To  complete  the  pictiire  we  consider  the  family  of  lines  through  0'. 

This  family  of  lines  includes,  among  others,  the  lines  representing  the 

six  wave  fronts.   Let  9   be  the  angle  that  the  normal  to  a  line  of  the 

family  makes  with  the  vertical  W-axis  through  0'.   In  Figure  10,  for 

example,  9     is  the  angle  that  the  slow  wave  front  makes  with  the  N-axis, 

Then,  it  is  easy  to  verify  that  (l)  the  pedal  curve  of  each  line  of  the 

family  consists  of  a  single  point  with  the  polar  coordinates  (ci.,9)  = 

(-  sin  0,  9)  and  (2)  the  locus  of  these  points  is  a  circle  of  diameter 
b 

V/b  whose  equation  is 

q  =  ^  sin  0.  (5.25) 

If  we  invert  this  circle  in  the  unit  circle  -  i.e.,  if  we  replace  q_ 
by  l/p,  we  find  that 

p  sin  9  =  ^  ,  (5.26) 

which  is  precisely  the  equation  of  the  line  £   of  equation  (5.22)  when^ 
as  is  the  case  here  a  =  0. 

C.   The  Graphical  Method 

Although  Huygens'  construction  has  greater  physical  appeal,  the 
graphical  method  outlined  at  the  end  of  subsection  A  has  much  in  its 
favor.   First,  it  is  more  accurate;  for  a  point  of  intersection  is 
more  easily  located  than  a  point  of  tangency.   Second,  as  the  angles 
of  reflection  and  refraction  and  corresponding  (reciprocal)  dis- 
tirrbance  speeds  are  determined  simultaneously,  it  yields  more  informa- 
tion.  Finally,  the  avoidance  of  the  cusped  figures  should  be  noted. 
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In  Figures  12(a)- (c)  we  show  polar  plots  for  r  =  l/2  (Figure 
12(a),  r  =  1  (Figure  12(b))  and  r  =  2  (Figure  12(c)).   In  all  cases,  a  =  0 
Q^   =  ii5  ^^^  i'-rr  =   0.   Figure  12(c)  also  contains  curves  associated 
with  the  parameter  values  r  =  2,  9     =  k'^     and  ^„  =  U5  ;  we  shall  have 
occasion  to  refer  to  these  later.   As  in  Figure  8,  all  angles  are 
measured  from  the  vertical  N-axis  and  are  considered  positive  or 
negative  according  as  they  are  measured  clockwise  or  counterclock- 
wise from  W. 

Whatever  the  value  of  r,   the  oval-like  curves  (labelled  f  ) 
correspond  to  the  fast  p.'  branch;  the  pair  of  solid  lines  (labelled 
A  ),  perpendicular  to  the  H_.-axis,  correspond  to  the  Alfven  p.  p- branch 
and  the  hyperbola-like  curves  (labelled  s^)  correspond  to  the  p  -slow 
branch.   The  pair  of  interrupted  lines,  perpendicular  to  the  H  -axis 


are  the  loci  of  points  satisfying  the  equation. 

pb  cos  ^jj  |cos  (0-0jj)|   =   [(1  +r)/r]^  (5.2?) 


These  lines  serve  as  asymptotes  of  the  s  -curves.   The  f ^ ,  A  and  s^- 
curves  are  all  symmetric  with  respect  to  the  origin  and  with  respect 
to  the  H  -axis.   It  follows  directly  from  equation  (5.1?)  and  the  cor- 
responding  equations  for  the  c.  ./b  that  when  ?/'  =  0  the  Alfven  branch 
intersects  the  slow  branch,  the  fast  branch  or  both  according  as  r 
exceeds  unity,  is  less  than  unity  or  equal  to  unity.   All  branches  have 
zero  curvature  on  the  H-p-axis  when  r  =[=  1.   When  r  =  1,  both  the  slow 
and  fast  branches  are  singular  on  the  H  -axis.   They  have  points  in 
common  at  p  =  1,  9=9     and  p  =  1,  9  =  9     +  n   and  at  each  of  these 
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points  are  tangent  to  a  pair  of  lines  that  make  an  angle  of  "^  = 
+  tan' "^(1/2)  with  the  H-axis.   [See  Figure  12(b)].   By  rotating  the 
curves  of  Figures  12(a),  (b),  or  (c )  about  the  H^  =  H-axis  we 
arrive  at  a  family  of  surfaces  in  p-space  known  collectively  as  the 
surface  of  phase  nomtials. 

The  curves  of  Figures  12(a)- (c)  are  typical  of  those  in  the 
parameter  classes  r  >  1,  r  =  1  and  r  <  1,  respectively.  Keeping 
r  and  Q„   fixed,  and  increasing  ^„,  we  find  from  the  equations  for 

n  n. 

the  c.  ./b  (cf.  equation  {h.k)-{h.'j))   that  for  each  &,   that  the  points 
on  the  f-,-ciurves  move  toward  the  origin  and  the  points  on  the  A  and 
So-cu3nres  move  away  from  the  origin.   The  general  appearance  of  these 
curves,  however,  remains  the  same  as  in  the  case  ■;//„  =  0  [see,  e.g., 
the  curves  of  Figure  11 (c)  associated  with  ^„  =  ho\    except  on  and 
in  the  neighborhood  of  the  H_-axis.   Here,  the  double  or  triple 
points  -  i.e.,  the  points  where  the  Alfven  wave  intersects  the  slow 
branch  (r  >  l),  the  fast  branch  (r  =  l)  or  both  (r  =  l)  -  split  up 
into  two  and  three  separate  points  respectively.   Moreover,  in  the 
case  r  =  1,  the  curvatures  of  the  slow  and  fast  branches  on  the  H 
axis  are  no  longer  singular.   Finally,  we  note  that  the  sole  effect 

of  varying  Q„  while  keeping  r  and  ^„  fixed  is  to  rotate  the  fp  A^ 

17 
and  s  -curves  rigidly  with  respect  to  the  fixed  axes  along  W  and  x_. 
3  -.     -^-u 


'Rather  than  draw  a  new  set  of  s,  A  and  f-ciirves  for  each  new  Q^,   it 
is  evidently  simpler  in  numerical  work  to  deal  with  one  fixed  set  of 
ctirves  and  to  rotate  the  W-axis  with  respect  to  this  fixed  set  in 
such  a  manner  that  H^.  makes  an  angle  0tt  with  N. 
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The  effect  of  varying  the  angular  parameters  ^„  and  9  may  be 
visualized  as  changing  the  orientation  of  the  surface  of  wave 
normals  with  respect  to  the  plane  of  incidence. 

Let  us  again  focus  on  the  problem  of  reflection  and  assume 
once  more  that  a  vanishes.   It  may  be  imagined  here  that  the  voa- 
disturbed  medium  in  [K  is   motionless  and  that  <xl'  is  a  rigid  surface. 
In  this  case^  as  we  pointed  out  earlier,  the  line  i   is  parallel  to 
the  N-axis;  its  distance  from  the  W-axis  is  evidently  the  reciprocal 
of  the  Alfven  number. 

By  varying  the  angle  of  incidence  and  counting  the  number  of 
intersections  of  £   with  the  lower  portions  of  the  f-,  ,   A  ,    s  -curves, 
one  arrives  at  the  following  conclusions:  (see  Figure  9a),  (l)  For 
sufficiently  small  angles  of  incidence  (and  correspondingly  large 
Alfven  numbers)  there  are  three  (real)  angles  of  reflection   ^--i^ 
0.p  and  9  (2)  To  each  fast  incident  wave  with  sufficiently 

small  1 9,  I  there  corresponds  three  such  angles  of  reflection.   How- 
ever, when  Q^    is  approximately  7O  the  fast  reflected  wave  is 

19   /  \ 
critically  reflected.     (3  J  When  the  incident  wave  is  of  the  slow 

or  Alfven  type  then  at  positive  and  negative  angles  of  incidence 

which  are  sufficiently  large  in  absolute  value,  the  fast  wave  is 

critically  reflected.  {h)   At  sufficiently  large  positive  angles  of 

incidence,  Alfven  or  slow  waves  can  produce  at  most  two  real  reflected 


__ 

Whether  or  not  the  associated  modes  of  propagation  are  excited  or  not 

depends  upon  the  boundary  conditions  and  on  the  orientation  of  the  mag- 
netic field. 

19 
■^See  footnote  I5. 
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waves  -  one  Alfven  and  one  slow  wave.   The  same  is  true  of  negative 
angles  of  incidence  which  are  moderate  in  absolute  value.   (5)  How- 
ever, at  negative  angles  of  incidence  that  become  successively 
larger  in  absolute  value,  first  the  Alfven  and  then  the  slow  wave 
is  critically  reflected. 

Statements  (l)-(5)  refer  specifically  to  Figure  9(a).   However, 
with  obvious  changes  they  apply  for  all  values  at  r,  ^„  and  9       the 
angles  ^„  =  +  9O  ,  +  I80  excepted.   Statement  (5)  is  no  longer 
valid  in  the  exceptional  cases;  but  (l)-(i|)  still  apply  with  minor 
changes.   The  reader  is  referred  to  reference  [23]  Figure  ^4  (b )  for 
numerical  results  and  in  addition  for  a  detailed  analysis  of  the 
problem  of  refraction. 

D.   General  Properties  of  the  Angles  of  Reflection 

The  angles  of  reflection  in  the  problem  discussed  above  enjoy 
several  general  properties  which  are  consequences  of  the  particular 
way  in  which  the  ratios  c.  ./b,  i,j  =  1,2,3,  depend  upon  9.,  r,  9 
and  ?//„.   We  wish  to  present  some  of  these  now. 

Let 

\j     =  ®ij  ^\'    'V  '^H^  ^]    ^'^   =1.2,3,  (5.28) 

express  the  functional  relationship  between  9. .,  9^    and  9   ,  f     and 
r.   Then  the  following  properties  can  be  shown  to  follow  from  Snell's 
laws  which,  in  the  present  case,  are 

sin  0j_  ^      sin  0^ 

i,j  =  1,2,3,  (5.29) 


c.  ./b        c./t 
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and  the  expressions  for  the  distrubance  speeds  (see  equations  (5.3)- 
(5.5)). 


©ij  [-'"J     =®ij  H  ^^^=1^3  A3 

®ij  l-\'  -^h]    =    -^ij,  \      (5.30) 

®ij  [-^i'  ^h1    =    -  ©ij  I^i'  -^H  ]    ^  S 

©,,  [-y    =-("-   Vi^^i>«  J  =  1.2.3.  A 

^"^  ^^  \-(Tt   +0.),    if  9.    <0; 


■6 


In  these  relations,  parameters  not  shown  are  considered  fixed  and 
unless  specifically  indicated  otherwise  the  subscripts  i  and  j  assume 
the  values  1,2,3.   The  angles  9,   .  appearing  in  equations  A^  and  A,^ 
are  what  we  shall  hereafter  refer  to  as  the  acute  angles  of  reflection. 
They  are  defined  as  follows: 


^i.1  =  \ 


e       if  0  >o    i  =  1,2,3, 
^^  ^^  (5.31) 


-(0.^  +:r),  if  0.^  <  0   j  =  1,2,3. 

Equation  A  expresses  the  fact  that  the  angles  of  reflection  are  in- 
dependent of  the  direction  of  the  component  of  H  which  is  perpendicular 
to  the  plane  of  incidence.   Equation  A  states  that  these  angles  are 
also  independent  of  the  sense  of  H^.   Together,  A  and  A  imply  that 
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the  reflected  angles  are  independent  of  the  sense  of  H.   Equation 
A  states  that  the  slow  and  fast  angles  of  reflection  are  invariant 
under  the  transformation  of  r  ->  r   whenever  the  incident  wave  is 
slow  or  fast.   Prom  A,  we  conclude  that  if  a  wave  incident  at  the 
angle  &.    produces  scattered  waves  at  the  angles  9.  .,   then  a  wave 
incident  at  the  angle  -6.   will  produce  scattered  waves  at  the  angles 
-9,   .  only  if  the  direction  of  H^.  is  reflected  in  the  W-axis.  specifi- 
cally,  only  if  9     is  replaced  hy  -^tt.   If  ®tt  is  not  so  replaced^  then 
according  to  A  the  scattered  waves  emerge  at  the  new  angles  9.  .   = 
-(ll'.  .  \j9 .;    ^u]  ;  which  differ^  in  general,  from  -Q .  ..      Only  when  9 
is  0  ,  +90  and  I80  do  we  find  that  9.  .   =  -0.  ..   In  other  words, 
negative  angles  of  incidence  are  not  physically  equivalent  to  positive 
angles  of  incidence  except  for  special  orientations  of  the  magnetic 
field. 

Equation  A^  expresses  the  reversibility  of  the  reflection 
process  in  the  following  sense.   If  9 .  .  is  the  acute  angle  of  re- 
flection  of  a  front  of  type  j  due  to  an  incident  front  of  type  i, 
then  a  front  of  type  j   incident  of  the  angle  -9,   .  will  produce  a 
reflected  front  of  type  i  which  emerges  at  the  angle  of  incidence  9, 
propagating  in  the  direction  opposite  to  that  of  the  original  in- 
cident front.   In  Aq,  it  is  assumed  that  the  9,   .  are  real.   These 

o'  ij 

relations  affirm  a  fact  to  be  expected  on  physical  grounds;  namely, 
that  the  fast  reflected  fronts  emerge  first,  followed  by  the  Alfven 
front,  and  then  the  slow  front . 
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The  analogues  of  A  -Aq  for  the  angles  of  refraction  in  the  case 
which  a'  vanishes  in  the  transmission  region  will  be  found  in 
reference  [23J  . 

E.   Explicit  Formulas  for  Alfven  Waves 
at  a  Contact  Discontinuity 

Suppose  the  surface  of  discontinuity  c?-^  is  a  hydroraagnetic  con- 

ri3i 

tact  discontinuity.   As  is  well  known ^  i-  -^  only  the  density  and 
entropy  may  "be  discontinuous  across  a  hydroraagnetic  contact  discon- 
tinuity, unless  HL,  the  normal  component  of  the  magnetic  field 
vanishes.   When  iL    vanishes,  a  discontinuity  of  u  across  9^   is  per- 
missihle;  however,  it  will  he  assumed  "below  that  u  is  continuous 
across  JO   in  this  exceptional  case  too.   Prom  the  Galilean  in- 
variance  of  the  jump  relations,  it  follows  that  u  may,  without  loss 
of  generality  "be  assumed  to  vanish  everywhere.   We  assume  below  that 
this  is  the  case. 

Suppose  now  that  an  Alfven  front  is  incident  at  the  angle  9^. 
We  wish  to  derive  the  acute  angle  of  reflection  0^^    (see  equation 

(5.31))  and  the  angle  of  refraction  0  „  at  which  the  scattered  Alfven 

?0 
fronts  emerge.    The  solution  of  this  problem  follows  almost  im- 
mediately from  Snell's  laws  (equation  (4.29))  and  equation  (5.5)  which 
implies  the  formulas  (0     ^  +  9^  ) 

ctn  022  =  "^^^  9^+2   tan  9^,  (5.32 


^°  When  H  is  not  in  the  plane  of  incidence,  slow  and  fast  fronts  will 
be  scattered  frome^.  Explicit  formulas  for  the  angles  corresponding 
to  these  waves  are  however  difficult  to  obtain. 


so  - 


ctn  e^^     =     ll^ctn  02  +  CH^"  D  tan  e^,  (5.33) 

where 

Tl  =  b/b-  =  (P'/P)''  (5.3^^) 

It  should  he  stressed  that  these  formulas  are  valid  for  all  values 
of  t//„  between  +  9O  ;  i.e.,  H  need  not  lie  in  the  plane  of  incidence.  It 
is  therefore  noteworthy  that  0pp  and  0p,  depend  only  on  9^.,  the  angle 
that  IL,  the  projection  of  H  on  the  plane  of  incidence,  makes  with  the 
normal.   This  interesting  fact  seems  to  have  been  overlooked  by  V. C.A. 
Ferraro  who  derived  another  somewhat  more  complicated  pair  of  equations. 

To  relate  Ferraro' s  equations  to  ours,  we  refer  to  Figure  7  and 
note  the  following  relations  between  0„  and  the  polar  angles  p  and  7 
employed  by  him: 

IL.  cos  0   =  H  sin  P 

H_  sin  0„  =  H  cos  P  cos  7, 
X.  n 

These,  in  turn,  imply  the  relation 

tan  0„  =  ctn  p  cos  7 
n. 

which,  combined  with  equations  (5.32)  and  (5.33)  yields  Ferraro's 

results  -  apart  from  obvious  notational  differences. 

In  Figure  I3  we  have  plotted  0pp  against  Qp  for  0„  =  0  , 

6     =  lj-5  and  0„  =  -ii5  .   These  are  typical  curves  in  the  parameter 
n.  11 

classes  0„  =  0,  0  >  0  and  0^.  <  0,  respectively. 

The  angle 

0*2  =  ctn"^  [  2  tan  0^^  ]  (5.35) 

introduced  in  Figure  13  is  the  maximum  positive  angle  of  reflection 

when  0„  is  positive  and  the  minimum  negative  angle  of  reflection 
xi 
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when  ^TT  is  negative.   For  a  given  ©„, 

02  =  ctn"^  [  -2  tan  9^^  (5.36) 

gives  the  critical  angle  of  reflection;  it  is  evidently  positive 

when  0„  is  negative  and  negative  when  9„  is  positive.   Wot  shown 
n  n 

in  Figiore  12  are  the  angles  of  reflection  associated  with  angles  of 

incidence  0     which  produce  anomalous  reflection  as  defined  in  the 

previous  section.   These^  to  he  sure,  roay  also  be  computed  from 

equation  (5.32)  and  (5.33).   Consider,  for  example,  the  curves  of 

Figure  13.   Here,  one  can  obtain  the  anomalous  values  of  9  _  by 

including  values  of  9     in  the  range  0^  <  0^  <  lt.5  when  9     =   -^5° 

and  by  including  values  of  0  in  the  range  9     >  9      >  -h"^     when 

0„  =  ^5  .   In  the  first  case  9         will  take  on  anomalous  values 
n  dd 

between  9O  and  135  while  in  the  second  case  9       will  take  on 
anomolous  values  between  -9O  and  -135  .   In  both  cases  the  re- 
flected fronts  will  be  found  in  the  same  quadrant  as  the  incident 
front.   The  reader  is  referred  to  reference  [22]  for  a  discussion  of 
equation  (5.33)  and  associated  graphs,  and  to  references  [7],  [22], 
[23J  for  a  derivation  of  the  coefficients  of  reflection  and  re- 
fraction. 


62  - 


APPENDIX  I 

Variation  of  the  Strengths  of  the  Fundamental  Waves 
Along  Their  Rays 

The  first  two  equations  in  the  recursive  matricial  system  (2.12)- 
(2.13)  are 

(M  -  p^c  I)  6R^°^  =  0  ,  (I.l) 

(M  -  p2  c  I)  5R  ^^^  =   (5R^°^)  ,  (1.2) 

where  6R  ^  '^  is  defined  in  equation  (2.I5).  Now  let  6R  ^  '^  be  any  one 

+  +    ±   ±    ±   ± 
of  the  six  solutions  g  R  ,   £»  R/,,   ^^  ^  ^  of  (l.l)  introduced  in 

Section  II  (see  equation  (2.28))  and  let  c  denote  the  corresponding 

disturbance  speed.   Since  M  is  symmetric  and  since  the  left  member 

of  (1.2)  has  the  same  matrix  as  the  left  member  of  (l.l),  it  follows 

that  (1.2)  has  a  solution  if  and  only  if  the  six  vector   P  (5  R^   )  is 

perpendicular  to  &  R    .   This  condition,  written  out  explicitly  with 

the  aid  of  equation  (2.I5)  reduces  after  some  manipulation  to  the 

orthogonality  condition  of  reference  [la,b] ;  namely, 

0  =  V  •  F  -.  p  5  J^^    •    (5  u^°^  .  V)u  -  i  (5P_^)'  (,  .  V)  ,2 
"-      ■*-        -       -       p       - 

(1.3) 

+  iM5H(°^)'v-  u-^.6h(°^.  (5H^°^.  v)u  . 

In  this  equation,  F  is  the  flux  vector  (see  equation  (2.5'))  associated 
with  the  given  fundamental  wave. 


63  - 


Note  that  5H^   ^  5u^  '^  and  6p^   are  linear  in  the  strength  e 
(cf.  equations  (2.2^1-)  -  (2.25))  of  the  wave.   We  may  therefore  cast 
equation  (l.3)  into  the  form 


V 


F  +  e  I  =  0 


(I.^) 


where  I  is  a  vector  which  is  independent  of  e.   Wow  we  have  seen 
(equation  (3.1^))  that  the  relation 

F  =  e2  V  ,  (1.5) 

V  denoting  the  ray  velocity;  obtains  in  each  of  the  six  fundamental 
waves.   Combining  this  relation  and  equation  (l.^),  we  find  that 

(1.6) 

But;  it  follows  from  simple  geometric  considerations  (see  equation 
(if.  6)  that 

d 


P  P 

v'Ve     +e      rV'V+ll      =     0. 
«,  i.  ,^  —I 


V  •    V     = 


dt 


log    (AU/U^) 


=  x(t)  J 


(1.7) 


Here,  A  is  the  expansion  ratio  defined  in  equation  (3.2l);  U  is  the 
normal  velocity  of  the  phase  front  at  time  t,  U  is  the  normal  velocity 
at  time  t  =  t  ana  x  =  x(t)  is  the  equation  of  the  ray.   Substituting 
this  expression  for  V  •  y  into  equation  (l.5)  and  noting  that  v  •  V  e 
is  the  derivative  of  e  along  x  =  x(t)  (i.e.,  that  v  •  V  e  =  de  /dt ) 
we  obtain  e  along  x  =  x(t)  as  a  solution  of  the  first  order  linear 
differential  equation 
2 


d  t 


I  + 


dt 


log  (AU/U^) 


X  =x(t)  J 


0. 


(1.8) 


Gh 


2 
The  expressions  for  €  catalogvied  in  Section  III  B  in  equations  (3.I7)- 

(3.20)  were  obtained  in  this  fashion. 

We  observe  finally  that  the  only  complex  quantity  in  equation  (I.8) 

is  £   (=  le  I  exp  (2i  (())  say).   It  follows  immediately  that  I 

d  t 

vanishes  along  the  ray. 
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FigureJ..   Cross- sectional  view  of  the  surface  of  nonnal  speeds  for  the 
parameter  value  r  =  I/2, 


s«^:(t) 
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s^7(t) 
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S^^lt) 


S^^(t) 
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Figure  2. 


Illustration  of  the  splitting  up  of  phase  front  element  bo/lt    ) 

f+  o' 

into  the  three  forward  propagating  elements  6  27.  (t),  i  =  s.  A, 

f  and  -ttiree  backward  propagating  elements  bS'    (t),  i  =  s.  A,  f. 


ee 


SX(t) 


>^(to) 


F  =  F(PJ 


ir(F) 


S^(to) 


/ 


/ 


H 


Figure  3.   Construction  of  the  fast  phase  front  J0>^(t)  evolving  out  of 
the  initial  front  Q/(t^)  in  the  case  where  the  unpertvirbed 
flow  is  spacially  uniform.   The  envelope  of  the  planes  T(F) 
determines  the  phase  front  at  time  t. 
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Jlgure  k.     The  Friedrichs'  ray  surface.  This  figure  depicts  the  phase 
fronts  that  emerged  from  a  point  source  at  0,  b   seconds 
earlier.  The  ray  surface  is  obtained  by  revolving  the  figure 
about  the  H-axis.  It  is  the  analogue  of  Presnel's  ray  surface 
in  optics.  The  undisturbed  medium  is  assumed  to  be  uniform. 
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Figure  '5.      I]J.ustration  of  the  fact  that 


and  (27  „  emerging  from  a 


f f 

spherical  source -of  (small)  diameter  2€  are  actually  distinct 

phase  fronts  or  order  2e   apart. 


dipole;^ 


Figure  6.  AlfVen  front  0  {O)   advancing  into  a  dipole  field. 

The  dipole  is  directed  along  the  z/d  axis;  d  is  the  fO^^I 
distance  of  90  {O)   from  the  dipole.   C  is  an  "earth" 
of  radius  l/^d.      Forward  propagating  fronts  are  shown 
in  the  upper  half  of  the  figure,  backward  propagating 
in  the  lower  half.   The  complete  picture  is  obtained 
by  adding  to  0  ~  (.1),  <v  7  (.2),  etc. ,  their  re- 
flections in  the  x/d  axis.  A  distinct  peak  forms 
at  "time"  t  =  .1+. 


-  7<J 


4z  =  N-axis 


H  cos/5  cos  y 
^X 


HsiniS 


Figure  7.  An  incident  wave  Iv  .    impinging  on  a  surface  of  discontinuity 

jy.  The  plane  I,  defined  by  the  normal  n.  to  the  incident  wave 
and  the  normal  N  is  the  plane  of  incidence.   Ht  is  the  pro- 
jection  of  H  on  the  plane  of  incidence.   The  angle  between  H 
and  H^  is  ^tt.  The  angle  (in  I)  between  IL.  and  N  is  &„. 
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Figure  8  (a).  The  disposition  of  the  phase  fronts  in  the  plane  of 
incidence. 
(b)  The  orientations  of  the  normals  to  the  phase  fronts  IV ^, 
tf/  .  and  ^.,  of  Figure  8(a). 
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(Pl3,^l3) 


Figure  9  (a)  Graphical  determinatiDn  of  the  angles  of  reflection.  AU.  angles  are  measured 

from  the  vertical  axis.   Clockwise  angles  are  regarded  positive,  counterclockwise  angles 

negative.  The  angles  of  reflection  9^.,   j   =  1,2,3  corresponding  to  an  angle  of  incidence 

0   of  a  front  of  type  i  are  obtained  hy  finding  the  intersections  of  the  line  i  with  the 
i' 


portions  of  s,  A  and  f  curves  that  lie  heneath  the  x-axis. 
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Figure  9,  (b)  Graphical  determination  of  the  arigles  of  reftaction  corresponding  to  an  angle 
of  incidence  0.  of  a  front  of  type  i.  The  figure  is  identical  with  the  upper  portion  of 
Figure  8(a).  However,  the  line  £'   that  determines  the  angles  of  refraction  is  closer  to 
or  further  from  the  vertical  axis  than  is  i  according  as  the  square  root  of  the  density- 
ratio  is  less  than  or  greater  than  unity. 
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Figure  11. 


The  projection  onto  the  plane  of  incidence  of  those  curves  on  the 
Friedrichs'  ray  surface  at  all  points  of  which  the  surface  normal 
is  parallel  to  the  plane  of  incidence.  The  parameters  r  and  9^ 
are  .5  and  0  respectively.   The  variable  parameter  is  '^„. 
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Figure  I3.   Plots  of  the  (acute)  angle  of  reflection  9       against  9     for 
the  parameter  values  9=0,   ajid  +  k^°. 
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